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Thermoacoustics  with  Multiple-Phase  Working  Fluids 


Abstract 

A  one  year  effort  to  improve  the  understanding  of  thermoacoustic  engines  and 
refrigerators  using  multiple-phase  working  fluids  resulted  in  a  more  complete.  The 
specific  system  considered  was  an  inert  gas-condensing  vapor  mixture  in  a  cylindrical 
pore  with  wet  walls.  It  was  found  that  vapor  diflusion  effects  in  the  mixture  are 
analogous  to  the  heat  diffusion  effects  in  the  thermoacoustics  of  inert  gases,  and  that 
these  effects  can  be  expressed  in  terms  of  the  thermoacoustic  viscous  function  used  in  the 
theory  of  sound  propagation  of  constant  cross-section  tubes.  This  theory  was  reduced  to 
a  concept  for  a  practical  device.  This  technical  report  in  contained,  for  the  most  part,  in 
two  papers  published  in  the  Journal  of  the  Acoustical  Society  of  America  which  are  a  part 
of  this  report  and  a  patent  also  included  here. 

1.0  Introduction 

Thermoacoustic  prime  movers  and  refrigerators  have  been  developed  over  several 
decades.’  These  devices  have  been  promoted  for  several  unique  advantages  over 
common  refrigeration  techniques,  particularly  in  low  incremental  construction  cost  and 
very  long  maintenance  cycles.  These  advantages  are  balanced  by  limitations  inherent  to 
acoustic  heat  engines  such  as  low  power  density,  low  efficiency,  and  complicated  heat 
exchangers.  Because  of  these  problems  acoustic  engines  have  not  made  any  substantial 
market  penetration  outside  of  orifice  pulse-tube  cryocoolers  on  spacecraft.^ 
Thermoacoustic  devices  use  acoustic  behavior  in  a  fluid  and  its  interaction  with  boundary 
surfaces  to  pump  heat.  A  refrigerator  will  move  heat  at  the  expense  of  acoustic  power, 
and  a  prime  mover  will  generator  power  from  a  heat  flow.  Most  machines  considered  or 
built  have  used  a  noble  gas,  usually  helium,  as  the  working  fluid.^  Occasionally  other 
simple  fluids  have  been  used,  including  mixtures  of  helium  and  another  noble  gas  or  a 
heavy  gas,'’  air,^  and  liquid  sodium.^  Machines  have  been  constructed  on  scales  ranging 
from  a  centimeter  to  tens  of  meters,  and  of  power  ranging  from  milliwatts^  up  to  tens  of 
kilowatts.®  In  every  case  the  efficiency,  defined  as  heat  pumping  power  to  input  power, 
and  power  density,  defined  as  heat  pumping  power  to  weight  or  volume,  have  been 
disappointing.  This  has  precluded  adoption  or  development  of  thermoacoustic  machines 
except  in  highly  specialized  applications.  The  theory  of  thermoacoustic  engines,  at  least 
in  the  linear  approximation,  has  been  well  established^’’®  ”  and  is  a  useful  starting  point 
for  engine  design. 

For  simple  fluids  the  linear  behavior  is  quite  well  understood.  Non-linear  effects, 
particularly  turbulent  dissipation  and  steady-state  streaming,  can  be  substantial  at  the  high 
acoustic  amplitudes  typical  of  thermoacoustic  engines,’^’®  and  are  generally  detrimental 
to  engine  performance.  While  these  nonlinear  effects  are  commonly  seen  and  discussed 
in  experimental  work,  they  are  difficult  to  predict  and  prevent.  The  predictions  reported 
here  are  based  on  linear  theories,  and  will  likely  overestimate  actual  performance.  Work 
extending  thermoacoustics  beyond  a  simple  working  fluid  has  developed  only  recently. 
Here  a  simple  fluid  is  defined  as  one  that  is  homogeneous  and  has  constant  composition. 


Experiments  with  heterogeneous  fluids  include  separation  of  noble  gas  mixtures  in  a 
thermoacoustic  device*'*  and  dehumidification  of  air  in  a  flow-through  thermoacoustic 
engine.*^  Analysis  of  these  devices  has  led  to  development  of  a  linear  theory  of 
thermoacoustics  in  heterogeneous  fluids*^  that  shows  a  large  change  in  heat  pumping 
behavior  when  part  of  the  working  fluid  is  evaporating  and  condensing  during  the 
acoustic  cycle.  These  theories  suggest  that  using  a  working  fluid  in  multiple  phases  may 
substantially  increase  the  power  density  and  efficiency  of  a  thermoacoustic  refrigerator 
over  a  similar  device  using  a  gaseous  working  fluid.  Such  an  enhancement  could  bring 
thermoacoustic  devices  to  the  range  of  power  density  and  cost  to  economically  compete 
with  current  technologies,  such  as  vapor-compression  cycles. 

In  order  to  test  these  ideas,  measurements  were  made  of  onset  temperature  in  an 
extremely  simple  thermoacoustic  prime  mover,  which  operates  in  atmospheric  air.*’  The 
condensable  fluid  used  is  water,  which  was  applied  to  the  stack  plates  with  an  atomizer. 
The  engine  is  a  modification  of  a  demonstration  device  developed  by  Reh-Lin  Chen.** 
Measurements  were  made  with  the  stack  at  various  positions  in  the  resonator,  and  with 
the  stack  either  dry  or  with  a  film  of  water.  These  experiments  showed  a  dramatic 
decrease,  by  about  a  factor  of  3,  in  the  temperature  gradient  required  for  the  onset  of 
thermoacoustic  oscillations  when  liquid  water  is  present  on  the  surface  of  the  stack. 

2.0  Technical  Accomplishments 

The  technical  results  of  this  study  are  contained  in  two  papers  and  one  patent  which  are 
attached  and  a  Master  Thesis  by  Dan  Brown  also  attached.  PhD  candidate  Tim  Simmons 
contributed  to  this  project  but  most  of  his  work  was  actually  supported  by  a  previous 
ONR  contract  so  his  dissertation  is  not  attached.  Copies  of  his  dissertation, 
“Experimental  Determination  of  Thermoacoustic  Stack  Properties”,  Aug,  2003  are 
available  upon  request. 
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Theory  of  inert  gas-condensing  vapor  thermoacoustics: 
Propagation  equation*^ 
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The  theory  of  acoustic  propagation  in  an  inert  gas-condensing  vapor  mixture  contained  in  a 
cylindrical  pore  with  wet  walls  and  an  imposed  temperature  gradient  is  developed.  It  is  shown  that 
the  vapor  diffusion  effects  in  the  mixture  are  analogous  to  the  heat  diffusion  effects  in  the 
thermoacoustics  of  inert  gases,  and  that  these  effects  occur  in  parallel  with  the  heat  diffusion  effects 
in  the  wet  system.  The  vapor  diffusion  effects  can  be  expressed  in  terms  of  the  thermoviscous 
function  F(X)  used  in  the  theory  of  sound  propagation  of  constant  cross-section  tubes.  As  such, 
these  results  can  be  extended  to  any  shape  parallel-walled  tube.  The  propagation  equations  predict 
that  the  temperature  gradient  required  for  onset  of  sound  amplification  in  a  wet-walled  prime  mover 
is  much  lower  than  the  corresponding  temperature  gradient  for  an  inert  gas  prime  mover.  The  results 
of  a  measurement  of  the  onset  temperature  of  a  simple  demonstration  prime  mover  in  air  with  a  dry 
stack  and  with  a  stack  wetted  with  water  provide  a  qualitative  verification  of  the  theory.  ©  2002 
Acoustical  Society  of  America,  [IX)I:  10.1121/1.1508113] 

PACS  numbers:  43.35.Ud  [MRS] 
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1.  INTRODUCTION 

Recent  worf:  has  developed  the  theory  of  sound  propa¬ 
gation  in  inert  gas-condensing  vapor  mixtures  in  narrow 
pores  with  wet  walls.  The  assumptions  made  about  the 
gas,  vapor,  and  pore  wall  are  straightforward,  namely:  the 
wavelength  and  characteristic  pore  dimension  must  be  much  ^ 
larger  than  the  mean  free  path  of  the  gas  and  vapor  mol¬ 
ecules,  steady  flow  of  the  gas  or  vapor  is  second  order  in 
acoustic  quantities  and  may  be  ignored  in  the  first-order 
analysis,  the  amplitude  of  the  acoustic  pressure  fluctuations 
is  Small,  and  end  effects  may  be  ignored.  It  is  assumed  that  a 
thin  layer  of  condensed  vapor  coats  the  solid  surface  of  the 
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pore  wall.  The  layer  is  assumed  thin  enough  to  neglect  its 
heat  capacity.  This  layer  is  the  site  for  evaporation  and  con¬ 
densation  of  vapor.  The  pores  are  assumed  to  be  on  the  order 
of  a  few  mass  diffusion  lengths  in  radius  so  that  condensa¬ 
tion  will  occur  preferentially  on  the  walls  rather  than  at  any 
condensation  centers  in  the  bulk  of  the  mixture.  Considering 
an  evaporation  and  condensation  process  implies  that  the  va¬ 
por  can  enter  the  thin  layer  but  the  inert  gas  cannot.  It  is  also 
assumed  that  there  are  no  other  sources  or  sinks  of  vapor  or 
inert  gas  within  the  pore,  that  the  heat  capacity  and  thermal 
conductivity  of  the  pore  walls  are  high  enough  to  ignore 
temperature  fluctuations  in  the  solid,  and  finally  that  the 
walls  can  be  considered  rigid. 

The  theory  demonstrates  that  the  terms  describing  the 
diffusion  of  the  vapor  are  analogous  in  form  to  die  terms 
describing  thermal  diffusion  in  the  mixture,  so  that  heat  and 
mass  transfer  terms  in  an  inert  gas-condensing  vapor  mixture 
in  a  thermoacoustic  engine  have  similar  forms.  If  mixtures 
could  be  found  vrith  heat  and  mass  transfer  terms  acting  to¬ 
gether,  then  the  heat  pumping  capacity  or  coefficient  of  per¬ 
formance  of  a  thermoacoustic  refrigerator  could  be  in¬ 
creased. 

The  thermodynamic  cycle  an  inert  gas-vapor  mixture 
undergoes  within  a  standing-wave  thermoacoustic  refrigera¬ 
tor  is  similar  to  that  of  a  pure  inert  gas,  with  the  addition  that 
during  the  compressive  part  of  the  acoustic  cycle  the  partial 
pressure  of  the  vapor  is  raised  above  the  saturation  pressure 
at  the  temperature  of  the  nearby  stack  wall  so  the  vapor 
condenses.  Conversely,  during  the  rarefaction  part  of  the 
acoustic  cycle  liquid  on  the  stack  evaporates.  In  this  way,  an 
inert  gas-vapor  working  fluid  will  transport  both  heat  and 
vapor  from  one  end  of  the  stack  to  the  other,  towards  the 
pressure  antinode.  If  a  temperature  gradient  is  imposed  such 
that  the  temperature  of  the  wall  increases  toward  the  antin¬ 
ode,  thermal  energy  will  be  transported  toward  the  antinode 
unless  the  temperature  of  the  wall  is  raised  to  the  point  that  it 
matches  or  exceeds  the  temperature  of  the  mixture.  This  con¬ 
dition  establishes  a  maximum  temperature  gradient  for  heat 
transfer  in  that  direction.  In  the  system  with  wet  walls,  in¬ 
creasing  the  wall  temperature  gradient  also  increases  the 
saturation  vapor  pressure  gradient  at  the  wall.  As  the  vapor 
pressure  gradient  at  the  wall  is  increased,  acoustic  vapor 
transport  will  occur  until  the  vapor  pressure  at  the  wall  is 
equal  to  the  partial  pressure  of  the  vapor  in  the  transported 
parcel.  This  establishes  a  different  critical  temperature  gradi¬ 
ent  for  the  acoustic  vapor  transport.  The  transport  of  heat 
from  the  mixture  to  the  wall  depends  solely  on  the  transverse 
temperature  gradients  in  the  pore.  The  transport  of  vapor  to 
the  wall  depends  solely  on  the  transverse  gradient  of  the 
partial  pressure  or  concentration  in  the  mixture.  The  vapor 
transport  also  carries  heat,  since  evaporating  and  condensing 
vapor  transports  entropy. 

In  this  paper  we  develop  the  acoustic  propagation  equa¬ 
tion  in  a  wet-walled  pore  with  an  imposed  temperature  gra¬ 
dient.  This  is  an  extension  of  the  dry  inert  gas  thermoacous- 
tic  propagation  formulation  known  as  Rott’s  equation.  This 
equation  will  be  analyzed  to  display  the  potential  of  inert 
gas-condensing  vapor  to  improve  thermoacoustic  prinie 
movers  and  refrigerators.  We  employ  the  notation  of  Amott 


et  ai^  for  the  thermoviscous  functions.  In  Sec.  R  the  hydro- 
dynamic  equations  describing  the  system  are  introduced  and 
Tijdeman’s^  “low  reduced  frequency”  approximation  is  ap¬ 
plied  to  this  system  of  equations.  The  gradient  in  ambient 
quantities  associated  with  the  ambient  temperature  gradient 
is  evaluated.  The  acoustic  boundary  conditions  are  described 
and  the  system  of  equations  is  then  solved  for  the  acoustic 
wave  equation.  In  Sec.  Ill  the  derived  form  of  the  wave 
equation  is  compared  to  the  form  found  in  Raspet  et  al}  for 
wet-walled  tubes  without  an  ambient  temperature  gradient 
and  to  the  form  for  dry  thermoacoustics."^’^  Section  IV  con¬ 
tains  a  discussion  of  our  formulation  with  respect  to  a  ther¬ 
moacoustic  prime  mover.  Section  V  contains  a  summary  of 
results  underlying  our  model  of  “wet”  thermoacoustics  and 
the  conclusions. 

II,  THEORY  OF  SOUND  PROPAGATION 
A.  Basic  equations 

The  model  presented  here  employs  the  assumptions  used 
by  Raspet  et  al}  in  the  calculation  of  sound  propagation  in 
wet  cylinders.  A  pore  consists  of  a  rigid  cylindrical  tube 
filled  with  an  inert  gas  and  saturated  vapor.  A  thin  layer  of 
liquid  coats  the  interior  of  the  tube.  It  is  assumed  that  the 
heat  capacity  of  the  tube  wall  is  high  enough  that  the  tem¬ 
perature  of  the  tube  wall  and  the  layer  of  liquid  do  not 
fluctuate.^  The  gas  component  carmot  penetrate  the  tube  wall 
or  the  liquid  layer,  and  the  gas-vapor  mixture  does  not  slip 
with  respect  to  the  liquid  layer.  The  gas  and  vapor  are  treated 
as  ideal  gases,  and  properties  of  the  mixture  are  computed 
according  to  the  theory  of  gas  mixtures.^  It  is  also  assumed 
that  terms  containing  ^7-,  the  thermal  diffusion  ratio,  and  kp , 
the  pressure  diffusion  ratio,  can  be  neglected.  The  effect  of 
kj  on  sound  propagation  in  wet  pores  has  been  shown  to  be 
small  for  this  system.^  The  effects  of  streaming  on  the  acous¬ 
tic  propagation  are  assumed  to  be  negligible. 

Equations  governing  sound  propagation  in  an  ideal  gas- 
vapor  mixture  are  presented.  The  Navier-Stokes  equation 
for  the  mixture  is 

D\  I  ii\ 

p  — =  -VP+/iV^v+l  )8+ jlV(V-v).  (1) 

The  continuity  equation  for  the  gas  can  be  written  in  terms  of 
the  gas  concentration  as 

P^+Vi=0,  (2) 

where  the  gas  concentration  is  C\—p\lp  and  the  density- 
weighted  species  difference  velocity  is  given  by  i 
“(PiP2^P)(vi““V2).  The  continuity  equation  for  the  vapor 
is 

^  +  V-(P2V2)=0,  (3) 

and  the  continuity  equation  for  the  mixture  is 

^  +  V-(pv)=0,  (4) 

where  the  mixture  density  is  given  by  p  =  Pi  +  P2 
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hydrodynamic  velocity  for  the  mixture  is  defined  as  v 
=  (piVi  +  P2V2)/p.  Only  two  independent  continuity  equa¬ 
tions  are  required  to  complete  the  system  of  equations.  The 
entropy  equation  for  the  mixture  is 

where  the  heat  flux  is  ^ven  by 


P={p0cl/y)(l+P*e-“^, 

(9) 

p=po(l+p*e"'"0. 

(10) 

C,=C^-bCfe"'"', 

(11) 

r=ro(i+r"e“'"0. 

(12) 

v=c,v*e"'"'. 

(13) 

•iSl- 


dT 


+  Mc 


P.c, 


the  chemical  potential  for  the  mixture  is 


III  fl2 
mi  m2 


i-/cvr,  (5b) 


(5c) 


where  fi'i  and  112  are  the  chemical  potential  per  molecule  of 
the  gas  and  vapor  respectively  and  tj'n.  is  the  viscous  stress 
tensor.*  The  ideal  gas  equation  of  state  is  used  for  both  the 
gas  and  vapor  components,  and  for  the  mixture 

P=nkT,  (6a) 

where  the  temperature  of  the  components  are  equal,  T=Ti 
=  T2,  the  total  mixture  pressure  is  the  sum  of  the  partial 
pressures,  P= ?/+  P2.  and  the  total  number  density  is  given 
jisn=ni+n2.In  terms  of  concentration  the  equation  of  state 
is 


P=nkT=pkT\ 


m,ffi2  /  m2 


The  diffusion  equation  for  the  mixture  is 

k 


i=-poZ>i2 


kr  icp 
VCi  +  ^VT+ 

^0 


VP 


(6b) 


(7) 


For  the  wet-wall  boundary  conditions  the  concentration  term 
dominates  the  diffusion  process  and  kj  and  kp  may  be  ne¬ 
glected.  This  differs  from  the  acoustically  driven  mass  diffu¬ 
sion  studies,^’’®  which  have  dry  pore  walls  and  the  principal 
contribution  to  diffusion  is  the  kp  term. 

The  ratio  of  specific  heats,  %  for  the  mixture  can  be 
related  to  the  ratios  of  the  two  component  gases  by 


ny 


A=v,-V2=c,A*e-'"',  .  (14) 

where  c,  denotes  the  ambient  speed  of  sound  in  the  mix¬ 
ture.  Symbols  with  a  superscripted  denote  small  di¬ 
mensionless  quantities,  while  a  subscripted  “0”  denotes 
ambient  quantities.  For  convenience,  we  define 

(15) 


v*=u*e^+v*er, 

and 

A*=U*e,+V*er. 

(5)  Define  the  normalized  axial  coordinate  as 
OfZ 

and  the  normalized  radial  coordinate  as 


where  R  is  the  tube  radius. 

(6)  Another  normalized  axial  variable  is  given  by 

z 


(16) 


(17) 


(18) 


(19) 


where  L  is  the  stack  length  and  will  be  used  to  describe 
ambient  gradients. 

The  components  of  the  Navier-Stokes  equation  are  un¬ 
changed  from  the  analysis  of  Raspet  et  al} 

{d^u^  1  du^  IcoRyd^u*] 


lU^ 


1  dP* 

y  H 

PqO)R^ 

dr) 


7)  dr] 


w 


y-1  7i-l  r2”l 

Please  see  the  list  of  symbols  for  a  complete  guide  to  vari- 
ables  and  symbols  used  in  the  equations. 

To  reduce  these  equations  to  a  useful  form,  we  follow 
the  analysis  of  Raspet  et  alP  For  the  sake  of  brevity  we  list 
the  required  steps  and  the  resulting  equations. 

(1)  Linearize  equations  assuming  that  steady  flows  are  neg¬ 
ligible  and  retaining  terms  in  dT^^ldz  and  dp^/dz^ 

(2)  Neglect  temperature-driven  diffusion  (Soret  effect)  and 
pressure-driven  diffusion  (Dufour  effect);  kj'=0f  kp^O. 

(3)  Assume  that  the  tubes  are  small  enough  that  we  may  use 
cylindrical  coordinates  with  no  azimuthal  dependence. 

(4)  Introduce  normalized  variables  with  time  depen¬ 
dence  \  ^ 


7  a>R\  dv*  u* 

(8) 

po(oR^\p,  3/^1 

[(77/  ^  V. 

(20) 


cAdP* 


y\o)Rl 

1  du^  u*  (cjRY^v^ 


1  du*  u*  lo)R\ 

\/t  3/577|\  c,  /  drj  rj 

The  continuity  equation  for  the  mixture  is 

i(— I 


/w/?\  du* 

1  d 

—  ___ 

V  dij 

u* 

R 

dpo 

Po 

L 

dZ 

(21) 


(22) 
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and  the  continuity  of  gas  in  terms  of  concentration  is 


P?P2fl  d  (<oR\dU* 


_| - - -  +M^ - - 

pQ  Po  ^  Po 


x/^i£-r  I  — 

To  °Ur\ro, 


Use  of  the  continuity  equations  and  the  diffusion  equation  in 
the  equation  of  state  yields 


Following  the  development  of  Landau  and  Lifshitz,*  Chap. 
VI,  and  using  the  ideal  gas  relation  for  the  thermal  expansion 
coefficient,  the  entropy  equation  is 


y—  1 

1 _ p* 

r 


K  l<aR\[^T*  \  dT* 

— -  - J  I  I  - 4*  “•  """" —  + 

PoCptoR^  \c,}[dij^  V  drj 

u*\R]dTf)  p\plu*\R\d 

To  [ij  rfZ  "pT 


771  ~W 


{P*-T*-p*} 


n®«2 

'*oPo 


a>R\  dU*  1  d 
Cs!  H  vdv 


^  ^.Po/>”l~'”2\  £  P1P2 

no  I  W|ffi2  j[L\\  dz[po\  Po  dz[  po 


The  axial  component  of  diffusion  is 


iei£iu*=. 

Po 


PoOia/PlPi/^o^^  ^fl 

Csjd^rj  dr]  ^  (  c,  / 


PoD,2/[1__^ 

wR‘-  \  [  9  <?!/  ^  \c,  j  d^' 


R  d  P?P2 

^LdZ~7r 


PoDi2/t/*[/?]^  d^  [pVz]  R^l<»R\  d  Ip^ip^ 

\  Po  L.  ^ .  Po  .  ^  Po  I  c,  /  d^  L  dz[  Po 


f/?l^</^fp?1  Ja)R\du*\R]  d\p° 
.*  _  _  Li  +  p2  _  _ _ Li 

l\  d^ipoj  \csjs^  l^^J^2:[po 


and  the  radial  component  of  diffusion  is 


0  0 

.PlP2_ 
l — T~^  ' 

Po 


^12  /  P1P2  [1  <?  ,  *, 


1  ^f/*[/?i  fp?pi 


Po  Po 


w/?\  dU*' 


<oR^\dv[L\dZ[po\l‘  '  ^ 

The  eight  equations  above  form  a  complete  set  for  the  vari¬ 
ables,  u*,  u*,  P*,  p*,  U*,  y*,  r*,  Cf .  We  note  that  the 
system  of  equations  can  be  solved  for  u*,  P*,  p*,  t/*, 
y*,  P*  without  using  Eq.  (23).  This  equation  is  retained 
here  for  convenience  in  future  calculations. 


vapor  in  the  mixture.  Since  the  total  pressure  in  the  resonator 
is  a  constant,  then  at  high  ambient  temperature  there  is  a 
higher  ratio  of  vapor  to  gas  in  the  mixture;  likewise,  at  lower 
ambient  temperature  there  is  a  lower  ratio  of  vapor  to  gas  in 
the  mixture. 

We  assume  that  the  pore  is  small  enough  so  that  the 
ambient  partial  pressure  is  a  function  only  of  the  axial  posi¬ 
tion.  This  is  equivalent  to  the  assumption  that  the  ambient 
gas  temperature  is  a  function  only  of  the  axial  position  as¬ 
sumed  in  inert  gas  thermoacoustics. 

Using  the  expressions  po==Pi  +  p2»  Pi^'ti'Wi,  and  p2 
—n2m2y  derivatives  with  respect  to  the  ambient  density  can 
be  written  as 


B.  Evaluation  of  ambient  gradients 

Derivatives  with  respect  to  the  ambient  density,  po,  and 
ambient  concentration,  p?/po,  must  be  evaluated.  In  dry 
thermoacoustics  the  ambient  density  is  inversely  propor¬ 
tional  to  the  ambient  temperature.  In  a  wet- walled  stack  with 
a  temperature  gradient,  the  temperature  dependence  of  the 
vapor  pressure  of  the  liquid  coating  the  stack  must  be  takeh 
into  account.  This  will  affect  the  ambient  ratio  of  gas  to 


driQ  dri2 

=mi— +  (m2-mi)  — 


driQ  dT^ 


dn\  dTo 


The  term,  driQ/dT^^  can  be  calculated  using  the  ideal  gas 
law,  yielding  dnQldT^=—nQlT(^  after  simplification.  The 
Clausius-Clapeyron**  equation  for  the  vapor  pressure  is 
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Ro\T^ 


=n2kTo, 


where  is  the  boiling  point  at  pressure  P^.  Evaluating 
for  dn\ldT(i  yields 

drA  Al  I  \ 


—  -  - - 1  .  (30) 

rfro  ToUoro  I 

We  now  introduce  the  “latent  heat  parameter,” 
(p=[(y-l)/r](l/Ro’Ro)r  as  defined  by  Marble.*^  Substitut¬ 
ing  into  Eq.  (28^  the  derivative  of  the  ambient  density,  after 
some  manipulation,  can  be  expressed  as 


dpo _ j_ 

dZ  ~  To  dZ  [  ^ 


«l\  P2 


y- 1  \  mil  Po 


Derivatives  of  the  ambient  gas  concentration,  P|/po>  can  be 
evaluated  in  a  similar  fashion 


A[£il=l^+  oIlL 

dZ  Pol  Po  dZ  '^^^dz{po. 


Mr  0^  o/p?  M  M 
z  (Pt+P2)-5z-Pi-^-PijZ 


mim2\  odf^odTo  dn^dTo] 
~~^[”^dTo  dZ  ”°dro  dZ\-  ^  ’ 

Making  use  of  the  derivatives  found  earlier  for  d«o  fdTo  and 
dn\ldTo,  our  expression  simplifies  to 

d(fi  dfp?l  y  p\p\no\dTo 


-dz=dz[Vor^y-^^ATodz- 

If  there  is  no  vapor  present,  then  diere  is  no  ambient  gradient 
in  the  gas  concentration.  Note  also  that  this  term,  like  the  one 
above  for  the  derivative  of  the  ambient  density,  is  scaled  by 
the  latent  heat  parameter. 


C.  Boundary  conditions 

The  boundary  conditions  at  the  wall  of  the  tube,  77=1, 
and  at  the  center,  77=0,  are  the  same  as  used  by  Raspet  et  al. 
These  boundary  conditions  must  be  cast  into  the  present  vari¬ 
ables.  The  boundary  conditions  at  the  tube  wall,  77=  1,  are  as 
follows. 

(1)  It  is  assumed  that  the  heat  capacity  of  the  tube  wall  is 
high  enough  that  the  temperature  of  the  tube  wall  and  the 
layer  of  condensate  do  not  fluctuate 

T*=0.  04) 

(2)  No  axial  slip  for  the  mixture  hydrodynamic  velocity 

M*  =0.  (35) 

(3)  The  gas  cannot  penetrate  the  tube  wall  or  condensate 
layer,  i;*=0,  which  can  be  written  as 

f,*+^V*=0.  (36) 

Po 

(4)  The  vapor  pressure  is  constant  at  the  wall  since  the  vapor 
pressure  is  only  a  function  of  temperature,  and  by  the 
first  boundary  condition  we  have  assumed  that  there  are 
no  temperature  fluctuations  at  the  wall.  This  implies  that 


^*  =  0  at  the  wall.  Use  of  the  continuity  conditions,  the 
definition  of  A,  and  the  equation  of  state  gives 


n?  l  d 


zftP*=--— (77V*). 
tiovdv 


The  boundary  conditions  on  the  aTcis  of  the  tube,  77=0, 
are  that  P*,  T*,  «*,  U*,  p* ,  and  p*  be  finite.  The  latter 
two  conditions  on  p*  and  p*  can  be  related  to  i>*  and  V*  by 
the  continuity  equations,  which  are  singular  unless  v  and  V 
are  zero  at  77=0.  Thus,  our  final  boundary  conditions  are 


(5)  Axis  of  tube  is  not  a  source  or  sink  of  vapor 

V*=0.  (38) 

(6)  Axis  of  tube  is  not  a  source  or  sink  of  mixture 

y*=0.  (39) 


D.  Solution 

A  formal  approach  to  the  approximation  and  interpreta¬ 
tion  of  the  complicated  equations  of  wet  thermoacoustics  is 
the  “low  reduced  frequency”  analysis  of  Tijdeman.®  This 
technique  has  been  applied  to  the  mass  transfer  problem  by 
Raspet  etaO 

The  principal  approximation  is  to  require  that  the  re¬ 
duced  firequency,  SI,  be  small 


(i)R  lirR 
Sl=  —  =  -— «^1. 

Cs  k 


(40) 


This  spans  a  wide  range  of  frequencies  for  small  tubes  and  in 
fact  includes  regions  where  the  Kirchhoff  wide-tube  approxi¬ 
mation  holds,  as  well  as  the  Rayleigh  narrow-tube/low- 
frequency  region.  The  physical  liieaning  of  this  approxima¬ 
tion  is  to  assume  that  the  tube  radius  is  much  smaller  than 
the  wavelength.  This  assures  that  there  wiU  be  no  cross 
modes  within  the  tube  and  that  radial  derivatives  are  much 
larger  than  axial  derivatives.  A  second  set  of  criteria  is  also 
used  to  retain  some  terms  that  would  otherwise  be  neglected 

u*>v*,  (41) 


and 

v*^V*.  (42) 

We  expect  the  first  condition  to  be  met  since  the  radial  ve¬ 
locity  must  be  zero  at  the  center  and  the  displacement  per 
cycle  is  at  most  on  the  order  of  one  radius.  The  second  con¬ 
dition  arises  since  v*  and  V*  are  the  same  order  of  magni¬ 
tude  at  the  wall  and  both  are  zero  at  the  center,  since  the  tube 
axis  is  not  a  source  or  sink.  The  above  criteria  are  the  same 
as  those  of  Raspet  et  al}  Generalization  to  include  an  ambi¬ 
ent  temperature  gradient  naturally  implies  a  length  scale, 
which  has  been  chosen  to  be  the  stack  length  L.  It  is  usually 
justifiable  to  assume  that  the  characteristic  pore  dimension, 
i.e.,  pore  radius,  is  much  smaller  than  the  stack  length 


(43) 


The  system  is  solvable  without  this  assumption,  but  this 
added  complexity  is  not  justified  in  this  study. 
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Certain  nondimensional  groupings  of  parameters  may  be 
identified; 

the  dimensionless  shear  wave  number 


l^_ 

>/2/? 

/  p. 

(44) 


(45) 


(46) 


the  dimensionless  thermal  wave  number 

Vnky  (o  Ipf^Cp  \{2R 

St  ’ 

and  the  dimensionless  diffusion  wave  number, 

Sd  ‘ 

These  wave  numbers  are  inversely  proportional  to  the  asso¬ 
ciated  penetration  depths. 

We  rewrite  Eqs.  (20)-(27)  using  Eqs.  (31)  and  (33)  for 
ambient  gradients  and  retain  the  lowest-order  terms  in  11  and 
RIL.  We  treat  terms  in  n*  and  V*  as  first-order  terms  in  ft, 
which  is  consistent  with  Eq.  (40)  and  Eq.  (43).  \Wth  the 
definitions  above  and  the  application  of  the  low  reduced  fre¬ 
quency  approximation,  the  axial  component  of  Navier— 
Stokes  is 


iu*= 


1  dP*  1 


r  K 


d^u*  1  Su* 
drj^  Tf  df) 


(47) 


and  the  radial  component  of  Navier-Stokes  is 

1  dP* 

--—=0. 

r  Srj 

The  axial  diffusion  equation  simplifies  to 
if/*  =  0, 

which  wiU  be  used  to  eliminate  U*  and  dU*l in  the  fol¬ 
lowing  equations.  The  radial  diffusion  equation  is 


(48) 

(49) 


1 

f'-D 


d7j\ri  d 7]^ 


y 

R 

1  dTo  du* 

<p  1 

r-1  rti 

L 

Tq  dZ  dr) 

(50) 


The  continuity  of  the  gas  in  terms  of  concentration  is 


0  0  -i  JJ 

iftCf  = — 5 ) 

•  pI  vSv 

0  0 

^  r  P1P2  '^0 

- T - T^  — 

v—  I 


r-i  Po  ”1 

the  continuity  equation  for  the  mixture  is 
du*  1  <?  , 

|•ftp*  =  ft•^+-T“(’7l>  ) 


lEl 

To  dZ' 


(51) 


To  dz\  r-l\  mijpo 
the  entropy  equation  is 
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(52) 


y—1  ft 


and  the  equation  of  state  is 

^0^0 

a{P*-T*-p*)=i-^imi-m2) 

”oPo 


—u*<p- 


d^T*  1  dT*' 

R 

u*  dTo 

Brf'  7}  dr) 

I 

To  dZ' 

(53) 


1  d 


y  ”0 

R 

1  dTo 

y-l  n\ 

I 

To  dZ 

(54) 


The  axial  and  radial  components  of  the  Navier-Stokes 
equation  are  unchanged  and  demonstrate  that  the  pressure  is 
not  a  function  of  radial  position.  The  axial  component  (47) 
shows  that  acoustic  velocity  is  generated  by  axial  gradients 
in  pressure,  with  some  modifications  to  account  for  shearing 
due  to  the  pore  walls.  The  terms  on  the  right-hand  side  of  the 
continuity  equation  (52)  represent  mass  entering  or  leaving  a 
region  due  to  acoustic  density  fluctuations  or  by  acoustic 
convection  of  a  density  gradient  The  entropy  equation  (53) 
describes  how  temperature  fluctuations  are  generated  by 
pressure  oscillations  and  by  fluid  oscillating  between  regions 
with  vaiying  ambient  temperature.  The  equation  of  state  (54) 
describes  how  the  density  fluctuations  are  generated  by  pres¬ 
sure  and  temperature  oscillations  as  well  as  diffusion  of  va¬ 
por  into  or  out  of  a  region.  Density  fluctuations  are  also 
generated  by  fluid  oscillating  back  and  forth  between  regions 
with  different  ambient  temperatures.  However,  if  the  two 
components  have  the  same  mass,  diffusion  does  not  contrib¬ 
ute  to  density  fluctuations  since  the  additional  term  describes 
density  changes  due  to  the  species  velocity  differing  from  the 
hydrodynamic  velocity.  Last,  the  radial  diffusion  equation 
(50)  describes  how  interspecies  diffusion  is  driven  by  gradi¬ 
ents  in  the  radial  diffusion  velocity  and  by  fluid  oscillating 
between  regions  of  different  ambient  temperatures. 

Note  that  the  Navier-Stokes,  entropy,  and  diffusion 
equations  are  decoupled.  The  solutions  will  be  in  terms  of 
separate  diffusion,  thermal,  and  viscous  wave  numbers.  This 
is  the  principal  result  of  the  low  reduced  frequency  approxi¬ 
mation.  This  shows  that  these  processes  are  independent,  to 
first  order,  and  only  coupled  in  the  equation  of  continuity  and 
equation  of  state  for  the  mixture. 

The  axial  Navier-Stokes  equation  is  solved  by  noting 
that  the  pressing  is  only  a  function  of  the  axial  coordinate 
and  by  applying  the  no-slip  boundary  condition  at  the  wall  to 
obtain  the  axial  hydrodynamic  velocity 


/  dP* 


1- 


idP* 


(55) 


We  have  introduced  here  the  definitions 


F,(r7)  =  l- 


where 
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F(\j)  is  th'ij  porous  media  dissipation  function  for  cylindri¬ 
cal  tubes  and  j=T,Dyfi/  The  recursion  relation  for  the 
Bessel  functions  has  been  used  to  evaluate  the  integral.  The 
evaluation  in  terms  of  generalized  dissipation  functions  is 
given  in  Ref.  4. 

Making  use  of  this  solution  and  the  boundary  condition 
on  vapor  pressure  at  the  wall,  the  radial  diffusion  equation  is 
solved  using  separation  of  variables  to  obtain  the  radial  dif¬ 
fusion  velocity 


11  <p  I  dTo\R]dP* 
ft  l-Nsc  ^  ~dz\L\dl 


«i-^scr-irodz 

MElJlkUhiHL] 
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where  the  Schmidt  number  is  fx! 

In  anticipation  of  later  results  it  is  convenient  to  calcu¬ 
late  the  fluctuating  component  of  the  gas  concentration,  C* 
by  making  use  of  the  solutions  for  m*  and  V*  above  to 
obtain 

0  0  "  r 

P1P2  ^0  r-  /  \\ 


1^1  dT^ 
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dP* 

ft  y-1  To  dZ 

I 

1^1  dFo 
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ft  y— 1  To  dZ 

I 

Multiply  through  by  rj dr}  and  integrate  from  0  to  1  using  the 
boundary  condition  on  y*  and  V*  at  7=0  to  eliminate  the 
constants  of  integration 


v*U=i= 


ift  d'^P*  F(X^)  ^  /ftp* 


-F{kr) 
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/  [P]  1  dPp  dP* 


[  1-Nf, 

Now,  employ  the  final  boundary  condition,  y*  +  (p2/po)^* 
=0  at  77=1,  using  the  above  expression  for  the  dimension¬ 
less  radial  diffusion  velocity  evaluated  at  ?7=1,  to  obtain  the 
following  form: 

d^p*  /  y-1  ”2  \ 

^F(X^)+  yP*[  1  -  ^F(\,.)+  4(1  -F(Xo)) 
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Converting  to  dimensional  form  in  terms  of  the  complex 
pressure  amplitude,  Pac~^o^*»  we  have 

d"^  P  AQ  dP  A  p  - 

-^+2a-^+k^PM:=0,  (63) 


Substituting  the  solution  for  u*  into  the  entropy  equation  and 
using  the  boundary  condition,  T*=0  at  17=!  yields  the  di¬ 
mensionless  temperature  fluctuations 

r*=— P*Fr(j;) 

y 

_J_  l_dP*[P]  1  <fFo[Fr(7;)-/Vi>rF«(^)' 

Cl  y  d^illTo  dz[  l-Np,  J’ 


where  the  Prandd  number  Npr=\^k^= fiCpI k. 

We  eliminate  p*  from  the  equation  of  state  using  Eq. 
(52)  and  substitute  the  solutions  for  T*,  V*,  and  «*  to  find 
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Equation  (63)  describes  the  acoustic  pressure  fluctuations  in 
a  wet-walled  tube  with  a  temperature  gradient.  We  note  that 
although  we  have  solved  this  for  cylindrical  pores,  parallel 
walled  pores  of  other  shapes  can  be  described  by  substitution 
of  the  appropriate  thermoviscous  dissipation  function, 
F(X;).'‘ 

For  purposes  of  discussion  and  numerical  integration,  it 
is  useful  to  rewrite  Eq.  (63)  as  two  first-order  differential 
equations  in  the  acoustic  pressure  and  volume  velocity.  Av¬ 
eraging  Eq.  (55),  the  axial  hydrodynamic  velocity,  over  the 
cross-sectional  area  of  the  pore  in  terms  of  the  dimensionless 
radial  coordinate  rj^rlR  and  ignoring  the  thin  condensate 
layer  gives 
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and 
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where  recursion  relations  for  the  Bessel  functions  have  been 
used  to  evaluate  the  integral.  Rearranging  to  solve  for  Ae 
acoustic  pressure  gradient  while  maldng  use  of  the  following 
expressions,  P*=Rac/^o»  Po~P(i^Vy>  ^—(ozlc^, 

and  realizing  that  the  acoustic  volumetric  velocity 

dP AC  _  ^voi  (67) 

dz  Anux  F{^ii)' 

where  Amix  cross-sectional  area  available  to 

the  gas/vapor  mixture.  Using  this  expression  in  Eq.  (63)  and 
using  the  expressions  above  for  P*,  P o .  a“d  ^  gives,  after 
rearrangement 

^  _  ^f(xjk2pj,c-2aU^,  (68) 

dz  lofpo 

where  a  and  are  defined  in  Eq.  (64)  and  Eq.  (65).  Equa¬ 
tions  (67)  and  (68)  arc  coupled  first-order  equations  for 
acoustic  pressure  fluctuations  and  acoustic  volumetric  veloc¬ 
ity  fluctuations. 


III.  LIMITING  CASES 
A.  Acoustics  of  wet  porous  media 

To  compare  to  recent  work  in  porous  media  theory,  we 
let  the  ambient  temperature  gradient,  dT^ldz,  go  to  zero  in 
Eqs.  (63)-(65).  This  makes  the  term  a=0,  and  our  acoustic 
wave  equation  becomes 

^^+k^PAC=0,  (69) 

dz 


with 
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F(X^) 


^f(x,)4(i- 
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This  agrees  with  Raspet  et  al}  subject  to  the  assumption  that 
the  thermal  diffusion  ratio,  kj,  is  zero. 


B.  Dry  thermoacoustics 

To  compare  to  dry  thermoacoustics,  we  let  the  number 
density  of  the  vapor,  n®-  8®  t®  2®™  (63)-(65).  The 

solution  for  the  wave  equation  then  becomes 
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Note  that  this  compares  with  the  solution  obtained  by  Amott 
et  aL^  To  compare  with  the  Rott  formulation,^  note  that  the 
relation  between  shear  wave  number,  and  viscous  pen¬ 
etration  depth,  ,  is  given  by  X^=  yJlR/ ,  and  the  rela¬ 
tion  between  thermal  wave  number,  Xy,  and  thermal  pen¬ 
etration  depth,  Sj<,  is  given  by  yJlR/Sfy  where  R  is  the 


ilc  iwirpi  thft  nore 


pore  perimeter.  The  relations  between  F(X^)  and  F{\t)  a®** 
the  Rott  functions,  and  f,^,  are  given  by,  F(X^)=1— 


and  F(Xr)= 1  •  The  complex  conjugate  arises  due  to  our 

use  of  e"'*"  notation,  rather  than  Rott’s  e'"'  notation. 


IV.  DISCUSSION 

Equations  (63)-(68)  clearly  demonstrate  that  the  ther¬ 
mal  diffusion  terms  for  inert  gas  thermoacoustics  are  aug¬ 
mented  by  additional  analogous  terms  in  mass  diffusion  in 
the  inert  gas-condensing  vapor  mixture  in  pores  with  wet 
walls.  Swift'^  discusses  the  contribution  of  each  term  in  the 
equations  of  inert  gas  thermoacoustics  corresponding  to  Eqs. 
(67)  and  (68).  The  term  corresponding  to  ~2aU\yt  is  iden¬ 
tified  as  the  temperature-gradient-dependent  complex  acous¬ 
tic  gain/attenuation.  For  standing-wave  devices,  the  imagi¬ 
nary  part  of  this  term  is  the  dominant  contribution  to  the 
gain;  in  a  pure  traveling  wave  device  it  is  the  real  part  that  is 
the  dorrunant  contribution.  Losses  arc  primarily  determined 
by  the  imaginary  part  of  k,  the  complex  square  root  of  Eq. 
(65).  For  the  inert  gas— vapor  mixtures  we  have  examined, 
X.D  and  \t  are  of  the  same  order  of  magnitude,  while  <p  can 
be  quite  large.  The  additional  gain  terms  are  larger  than  the 
additional  loss  terms  so  that  lower  onset  temperatures  and 
higher  woric  can  be  expected  from  a  wet  prime  mover. 

In  order  to  test  these  ideas,  measurements  were  made  of 
onset  temperature  in  an  extremely  simple  thermoacoustic 
prime  mover,  which  operates  in  atmospheric  air.  The  gas 
properties  needed  to  compare  the  dry  and  wet  gain  and  loss 
terms  have  been  estimated  using  the  methods  of  Ref.  1.  The 
condensable  fluid  used  is  water,  which  was  appHed  to  the 
stack  plates  with  an  atomizer.  The  engine  is  a  modifica^on  of 
a  demonstration  device  developed  by  Reh-Lin  Chen.  The 
resonator  section  is  a  lOO-ml  Pyrex  graduated  cylinder  (di¬ 
ameter  26  mm,  length  240  mm)  which  forms  a  closed-open 
resonator  exposed  to  the  atmosphere.  The  stack  is  a  1-in.- 
long  cylindrical  section  of  Celcor  square-celled  ceramic  hon¬ 
eycomb  with  200  cells  per  square  inch.  The  lower  end  is 
woven  with  about  250  mm  of  30-AWG  nichrome  wire  which 
can  be  electrically  heated.  The  upper,  cold,  end  of  the  stack 
has  no  heat  exchanger  at  all,  and  dissipates  heat  into  the 
nearby  glass  wall,  convection  of  the  air  above  it,  and  its  own 
heat  capacity. 

Beginning  with  the  entire  system  at  ambient  tempera¬ 
ture,  electrical  power  is  applied  to  the  heating  wire  over 


J.  Acoust.  Soc.  Am..  Vol.  112,  No.  4.  October  2002 


Raspet  et  aL:  Thermoacoustic  propagation  in  condensing  mixtures  1421 


20-30  s  t©  produce  an  increasing  temperature  gradient 
across  the  stack,  measured  with  thermocouples  cemented  to 
the  stack  at  each  end.  The  temperatures  are  recorded  when 
sound  is  detected  radiating  from  the  open  end  of  the  resona¬ 
tor.  Measurements  were  made  with  the  stack  at  various  po¬ 
sitions  in  the  resonator,  and  with  the  stack  either  dry  or  with 
a  film  of  water.  The  stack  was  located  near  its  optimum 
location  for  low  onset  temperature.  Since  there  is  no  cold 
heat  exchanger,  these  experiments  were  performed  with 
rapid  heating.  With  only  two  thermocouples,  we  could  not 
measure  the  actual  temperature  distribution  in  the  stack.  For 
this  reason  careful  comparison  of  data  and  theory  is  not  jus¬ 
tified.  These  experiments  showed  a  dramatic  decrease  in  the 
temperature  difference  required  for  the  onset  of  thermoa¬ 
coustic  oscillations  when  liquid  water  was  present  on  the 
surface  of  the  stack.  With  a  wet  stack,  the  onset  temperature 
difference  was  about  80  K,  with  a  dry  stack  about  280  K.  At 
the  average  temperature  of  the  wet  stack  (330  K),  n2//ii 
=0.22,  <^=4.2,  7=1.388.  This  gives  a  coefficient  of  about 
3.2  for  the  diffusion-related  gain  term  in  Eq.  (64).  The  coef¬ 
ficient  of  the  diffusion  loss  term  is  0.22  in  Eq.  (65),  while  the 
coefficient  of  the  thermal  loss  term  is  about  0.29.  The  vis¬ 
cous  losses  are  usually  larger  than  the  thermal  losses,  and  the 
viscosity  of  dry  air  at  430  K  is  larger  than  the  viscosity  of 
saturated  air  at  330  K  (2.4X  10”®  vs  1.8 10"®  kg/ins).  If  the 
viscous,  thermal,  and  diffusion  lengths  are  approximated  as 
equal,  then  the  ratio  of  the  critical  temperature  gradient  with 
and  without  diffusion  should  be  on  the  order  of  4.  This 
agrees  reasonably  well  with  the  measured  ratio  of  3.5.  The 
measurements  provide  qualitative  confirmation  that  the  gains 
predicted  by  the  theory  can  be  realized  in  a  real  system. 
Measurements  with  better  instrumentation  may  provide 
quantitative  confirmation  in  the  future. 

These  results  are  promising  for  improving  the  heat  flow 
and  efficiency  of  thermoacoustic  refrigeration,  since  the  mass 
diffusion  effects  act  in  parallel  with  the  heat  diffusion  effects 
and  can  augment  the  heat  transport.  To  analyze  such  a  system 
requires  expressions  for  the  enthalpy,  heat,  and  work  flows 
within  the  stack.  The  derivation  of  these  expressions  and  the 
analysis  of  “wet”  thermoacoustic  refrigeration  are  presented 
in  the  following  paper. 

V.  CONCLUSIONS 

The  basic  propagation  equations  and  equations  for  first- 
order  acoustic  variables  have  been  derived  for  an  inert  gas- 
condensing  vapor  mixture  in  wet-walled  tubes  with  imposed 


temperature  gradients.  These  equations  demonstrate  that  the 
addition  of  wet  walls  and  condensing  and  evaporating  vapor 
increases  the  acoustic  gain  produced  for  a  given  temperature 
difference  across  the  stack.  A  simple  measurement  demon¬ 
strates  that  the  predicted  gains  can  be  realized  in  physical 
systems.  The  derived  equations  also  form  the  basis  for  the 
development  of  woik  and  heat  flow  equations  and  therefore 
the  analysis  of  “wet”  thermoacoustic  refrigeration. 
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The  preceding  paper  [J.  Acoust.  Soc.  Am.  112, 1414-1422  (2002)]  derives  the  propagation  equation 
for  sound  in  an  inert  gas-condensing  vapor  mixture  in  a  wet-walled  pore  with  an  imposed 
temperature  gradient.  In  this  paper  the  mass,  enthalpy,  heat,  and  woric  transport  equations  necessary 
to  describe  the  steady-state  operation  of  a  wet-waUed  thermoacoustic  refrigerator  are  derived  and 
presented  in  a  form  suitable  for  numerical  evaluation.  The  requirement  that  the  refrigerator  operate 
in  the  steady  state  imposes  zero  mass  flux  for  each  species  through  a  cross  section.  This  in  turn  leads 
to  the  evaluation  of  the  mass  flux  of  vapor  in  the  system.  The  vapor  transport  and  heat  transport  ^ 
shown  to  work  in  parallel  to  produce  additional  cooling  power  in  the  wet  refrigerator.  An  idealized 
calculation  of  the  coefficient  of  performance  (COP)  of  a  wet-walled  thermoacoustic  refrigerator  is 
derived  and  evaluated  for  a  refrigeration  system.  The  results  of  this  calculation  indicate  that  the 
wet-walled  system  can  improve  the  performance  of  thermoacoustic  refrigerators.  Several 
experimental  and  practical  questions  and  problems  that  must  be  addressed  before  a  practical  device 
can  be  designed  and  tested  are  described.  ©  2002  Acoustical  Society  of  America, 
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pation  functions  in  Rott’s  notation.^  The  use  of  Amott’s  no¬ 
tation  is  consistent  with  our  earlier  papers  on  sound  propa¬ 
gation  in  wet  pores.^”^ 

In  a  “wet”  thermoacoustic  refrigerator,  energy  is  trans¬ 
ported  along  the  temperature  gradient  in  the  stack  both  in  the 
form  of  heat  and  as  energy  carried  by  the  condensable  spe¬ 
cies.  The  vapor  condenses  and  evaporates  on  a  thin  layer  of 
liquid  on  the  wall  of  the  pore.  Proper  phasing  between  the 
evaporation/condensation  of  vapor  from  this  film  and  the 
axial  acoustic  velocity  will  result  in  net  pumping  of  vapor 
from  one  end  of  the  stack  to  the  other.  In  its  simplest  form, 
we  suppose  that  the  condensable  fluid  returns  by  flow  in  tlie 
thin  liquid  layer.  This  represents  a  loss  mechanism  analogous 
to  heat  conduction  in  the  stack  and  its  contribution  must  be 
evaluated. 

DeltaE  is  a  well-established  numerical  code  for  model¬ 
ing  the  performance  of  thermoacoustic  devices.*  The  trans¬ 
port  equations  will  be  evaluated  in  forms  suitable  for  use  in 
DeltaE.  DeltaE  uses  the  conservation  of  the  total  time- 
averaged,  second-order  steady-state  enthalpy  flow  within  the 
stack  to  calculate  the  temperature  ^dient  and  temperature 
profile  through  the  stack.  The  enthalpy  is  conserved  because 
the  stack  and  refrigerator  are  assumed  to  be  insulated  and 
that  heat  exchange  with  the  environment  only  occurs  at  the 
heat  exchangers  at  each  end  of  the  stack.  This  assumption  is 
also  valid  for  the  wet  system. 

In  Sec.  n,  the  first-order  differential  equations  for  the 
volumetric  velocity  and  acoustic  pressure  are  presented  and 
their  use  in  numerical  calculations  discussed.  Section  III  de¬ 
velops  the  expression  for  the  mass  flux  in  the  stack  in  steady- 
state  operation.  Section  IV  develops  the  equations  for  the 
time-averaged  enthalpy  flux,  aind  the  heat  and  work  fluxes  in 
the  stack.  From  the  expression  for  the  total  enthalpy  flow,  the 


I.  INTRODUCTION 

In  a  companion  paper  [J.  Acoust,  Soc.  Am.  112,  1414- 
1422  (2002)],  hereafter  referred  to  as  Paper  I,  the  acoustic 
propagation  equation  was  derived  for  the  pressure  wave  in  an 
inert  gas-condensing  vapor  mixture  in  a  wet-walled  pore 
with  an  applied  temperature  gradient.^  The  expressions  for 
other  first-order  acoustic  variables  were  also  derived,  such  as 
the  acoustic  temperature  fluctuations,  acoustic  velocity  fluc¬ 
tuations,  and  acoustic  concentration  fluctuations.  A  complete 
analysis  of  the  potential  of  the  inert  gas-condensing  vapor 
system  to  thermoacoustic  refrigeration  requires  expressions 
for  the  total  axial  enthalpy  flux,  heat,  and  work  fluxes,  as 
well  as  the  total  mass  flux.  With  these  expressions  the  coef¬ 
ficient  of  performance  can  be  evaluated  and  compared  to  the 
analogous  expressions  for  inert  gas  thermoacoustics.  In  ad¬ 
dition,  the  system  of  equations  for  numerical  integration  of 
realistic  stacks  can  be  developed.  In  the  derivations  we  use 
Swift’s  review  article^  as  a  guide,  but  employ  the  notation  of 
Amott  et  al,^  Cylindrical  geometry  is  chosen  for  the  pores 
within  the  stack  for  simplicity.  The  resulting  equations  may 
be  adapted  to  other  pore  shapes  by  substituting  different 
functional  forms  of  the  porous  media  dissipation  functions. 
The  work  of  Amott  et  al,^  contains  expressions  for  the  po¬ 
rous  media  dissipation  functions  for  various  pore  geometries. 
Also,  Swift’s  thermoacoustics  textbook  contains  these  dissi- 


•^Portions  of  this  work  have  been  presented  in  William  V.  Slaton  and  Rich¬ 
ard  Raspet,  “Wet- walled  Thermoacoustics,”  First  International  Workshop 
on  Thermoacoustics,  ’s-Hertogenbosch,  The  Netherlands,  22-25  April 
(2001);  William  V.  Slaton  and  Richard  Raspet,  “Wet-walled  thermoacous- 
’  tics,”  J.  Acoust.  Soc.  Am.  110,  2677  (2001). 
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heat  and  mass  pumping  by  acoustic  processes  and  loss 
mechanisms  may  be  identified.  Section  V  describes  the  pro¬ 
cess  occurring  in  a  wet  thcrmoacoustic  refrigerator  stack  and 
identifies  a  useful  regime  of  operation.  The  coefficient  of 
performance  of  an  idealized  inert  gas-vapor  thermoacoustic 
re&igerator  is  developed  in  Sec.  VI.  This  expression  for  the 
coefficient  of  performance  is  then  analyzed  to  determine  the 
potential  of  wet-walled  system  to  improve  the  practicality  of 
thermoacoustic  refirigeration.  Section  VII  contains  the  con¬ 
clusions  of  this  research. 


II.  PROPAGATION  EQUATIONS 

DeltaE  numaically  integrates  the  two  coupled  first- 
order  differential  equations  in  acoustic  pressure  and  volumet¬ 
ric  velocity  to  determine  the  acoustic  propagation  through 
the  stack.  From  Ref.  1  we  have 


dP^c  »fe>Po  ^vol  (1) 

dz 

where  Ambc  is  folsl  stack  cross-sectional  area  available  to 
the  gas-vapor  mixture,  l/yoi  is  the  volume  velocity  in  the 
stack  governed  by 


dUyfH  ^ 

dz 


lOfpo 

where  a  and  are  defined  as 


F(X„)k2PAC-2«t^vol. 


(x~  X : 


1  dTo 

fF(Xr)/F(X^)-l 

To  dz  1 

y  F(Xp)/F(X^)-l 


y-1 


l-Ni 


Sc 


and 


<a^  y 


F(X  J 


(2) 


1  -  - — F(Xr)+ *4(1  ~F(X£))) 

r  «i 


(3) 


Equations  (1)  and  (2)  are  coupled  first-order  equations  for 
acoustic  pressure  fluctuations  and  acoustic  volumetric  veloc¬ 
ity  fluctuations.  In  steady  state,  knowledge  of  the  value  of 
the  fluctuating  variables  and  the  total  enthalpy  within  the 
stack  at  one  location  allows  calculation  of  the  fluctuating 
variables  at  another  location  within  the  stack.  It  can  be  seen 
that  thermal  and  mass  transfer  processes  have  similar  func¬ 
tional  forms  in  Eqs.  (3)  and  (4).  The  behavior  of  the  ambient 
temperature  gradient  is  evaluated  from  the  enthalpy  equation 
to  be  derived  in  Sec.  IV. 


111.  MASS  FLUX 

Within  a  closed  resonator  in  steady  state,  the  second- 
order  time-averaged  mass  flux  of  the  inert  gas  must  be  zero. 
Acoustic  pumping  of  inert  gas  in  the  axial  direction  must  be 
balanced  by  an  axial  return  mean  flow  of  gas  to  equalize  the 
mean  pressure  within  the  resonator.  Thus,  we  may  write 

j  (piy+p°ys,i)idA^=0,  (5) 


where  is  the  steady-state  flow  velocity  of  inert  gas.  v  is 
the  first-order  acoustic  fluctuation  velocity.  In  Ref.  I  it  is 
shown  that  there  is  no  axial  acoustic  difference  velocity  term 
in  the  low  reduced  firequency  solution.  Therefore,  the  acous¬ 
tic  velocity  of  each  species  is  the  same  as  the  mixture  acous¬ 
tic  velocity.  We  have  proposed  an  acoustic  propagation 
model  in  which  there  is  a  net  mass  flux  of  vapor  due  to 
acoustic  processes.  An  analogous  exjnession  may  be  written 
for  the  vapor 


J  ^vzpoc  ^liquid  >  (6) 

where  m^quid  is  the  mass  flow  rate  in  the  condensed  liquid 
layer  and  is  the  steady-state  flow  velocity  of  the  vapor. 
For  steady-state  operation  the  vapor  transported  in  one  direc¬ 
tion  in  the  mixture  must  return  as  liquid  flow  in  the  other 
diiection.  The  ambient  steady  flow  of  inert  gas  and  vapor  is 
also  related  to  the  ambient  concentration  of  the  inert  gas 
within  the  stack.  This  is  illustrated  by  recalling  the  expres¬ 
sion  for  the  ambient-density-weighted  species  difference  ve¬ 
locity 

0  0 

io=-poD,2VC?=-^(v,.,-v,2),  (7) 

Po 


and  that  the  acoustic  concentration  fluctuation  can  be  ex¬ 
pressed  by 


c:f= 


p?p5(pf-p|) 

- ^ - 


(8) 


Equations  (5)-(8)  may  be  combined  and  rearranged  to  form 
two  useful  expressions.  The  iSrst  is 


Pi  P? 

^12^  ^l)z^^inix~  ““  T  ^ vapor”"  ^liquid  ► 

Po  Po 


(9) 


and  the  second  is 


J  ^vapor  ^liquid*  (1®) 

where  v^^CpiV^j+p^v^^/po  is  the  ambient  hydrodynamic 

steady  flow  of  the  mixture. 

Using  expressions  for  C®,  Cf ,  and  from  Ref.  1  and 
expressing  the  time  average  of  the  product  of  the  acoustic 
quantities  as  yRe[w*Cf],  we  find 

^0  y  ^iiiix^l2P2  ^^0  /11^ 

where  ^=[(7“*  IVrK^/^oJ’o)- 

Substitution  of  w*  and  Cf  yield  integrals  of  F|(  rj)  or  of 
products  F/(77)Fy  (77).  Reference  4  shows  that  these  inte¬ 
grals  are  given  as 

f  F,(77)F,(,7)’7‘/’7=5Re[F(X,)],  02a) 

Jo 
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I. 


1  F(\j)Ni.+F(\,) 


and 


rFi(v)vJv=jF(k!), 

Jo 

where  A^y=Xf/\y. 

Using  these  integrals,  we  find 


(12b) 


(12c) 


1  0  1 
1  P2  no  Y  1 


2po^Y-iroC^ 


Re[PAcf/voiH,(Xp)] 


-<P 


Cppo  dip 

Uy<A 

<oA^  dz 

F(X^) 

Im[E2(XD)] 


np  y  AmhPiaPz  dT^ 

"^^lY-l  To  dz’ 


where  we  introduce  the  notation 


B,(xp= 


f(X,)/F(X^)-l 

l+AT,- 


(13) 


(14) 


integral  over  the  whole  stack  can  be  broken  up  into  integrals 
over  the  solid,  liquid,  and  open  areas  of  the  stack.  We  retain 
second-order  terms  in  acoustically  fluctuating  variables  and 
steady-state  terms  in  the  time  average,  and  ignore  acoustic 
terms  in  the  solid  and  liquid  layer  and  all  viscous  heating 
terms,  v-<r'.  Streaming  velocities  are  assumed  to  be  at  least 
first  order  in  acoustic  quantities.  Thus,  the  axial  component 
of  the  total  enthalpy  flux  is 

^=f  ('^Lx('^P+Po^s)+Poyw+<lo)zdA^ 


where  the  overbars  denote  time  averaging,  and  the  ambient 
loss  mechanisms  within  the  gas  mixture  are  contained  within 


qo= 


P.c, 


io-KVTo. 


(18) 


and 

»2(X;)  . 


(15) 


withy  taken  as  T referring  to  thermal  effects  or  D  referring  to 
mass  diffusion  effects  and  Nj  referring  either  to  the  Prandtl 
number  for  thermal  effects  or  the  Schmidt  number  for  mass 
diffusion  effects. 

The  acoustically  driven  terms  can  be  recognized  as  very 
similar  in  form  to  the  heat  transport  equations  for  inert  gas 
thermoacoustics:  there  is  an  acoustically  driven  term  inde¬ 
pendent  of  temperature  gradient  and  another  acoustically 
driven  term  proportional  to  temperature  gradient.  Determina¬ 
tion  of  the  direction  of  these  terms  is  most  easily  accom¬ 
plished  by  comparison  with  the  well-understood  enthalpy 
flux  terms  of  inert  gas  thermoacoustics,  so  we  defer  the  dis¬ 
cussion  of  the  direction  of  mass  flux  until  the  enthalpy  and 
heat  flux  terms  have  been  derived  and  discussed. 


The  Soret  effect  coefficient,  and  the  Dufour  effect  coef¬ 
ficient,  Jt/»,  have  been  ignored  to  conform  with  earlier  deri¬ 
vations.  fic  is  the  chemical  potential  of  the  mixture  and  is 
defined  as  where  and  are 

the  chemical  potentials  per  molecule.  Using  Eq.  (7)  for  the 
ambient-density-weighted  species  difference  velocity,  Eq. 
(18)  may  be  rewritten  as 

qo=-(H'?-w5)poO,2Vc?-Avro,  (i9) 


where  Wj  is  the  ambient  enthalpy  per  unit  mass  of  species  y. 

Equation  (10)  can  be  used  to  express  the  enthalpy  flux 
for  the  mixture  in  terms  of  or  myapor 


^  (^mix  ^ liquid) ^vapor*^  ^  | 


^  .  dC?  dTo\ 


dz 


IV.  ENTHALPY  BALANCE 

The  enthalpy  flow  expression  is  derived  from  the  energy 
density  equation  given  by  Landau  and  Lifshitz^ 

v-a'  +  qo  . 

(16) 

e  is  the  internal  energy  per  unit  mass,  w  is  the  enthalpy  per 
unit  mass,  v*cr'  is  the  energy  flux  due  to  internal  friction, 
and  qo  is  the  heat  flux.  Use  of  this  basic  expression  insures 
that  we  have  included  all  relevant  terms  in  our  derivation  at 
the  cost  of  initial  complexity.  The  terms  in  the  divergence 
represent  the  energy  flux  within  the  stack.  The  axial  compo¬ 
nent  of  the  enthalpy  flux  within  the  stack  is  of  interest  be¬ 
cause  it  is  a  conserved  quantity  in  steady-state  operation.  The 
time-averaged  axial  enthalpy  flow  is  found  by  averaging  Eq. 
(16)  over  the  whole  stack  and  taking  the  time  average.  The 


JaiQ  I  «i0 

^soUd  ^‘^solid  J  ^liquid  liquid*  (^^) 

Terms  in  Eq.  (20)  may  be  identified:  PqUW  represents  the 
enthalpy  flow  down  the  stack  due  to  acoustics.  The  term 
containing  Wj  — represents  the  enthalpy  flow  due  to  the 
ambient  diffusion  of  the  inert  gas  and  vapor.  The  term 
(»*'2iix“*^uquid)'”vapor  represents  the  enthalpy  flow  due  to  the 
net  transport  of  the  gas  mixture  and  the  return  flow  in  the 
liquid  layer.  The  diffusion  loss  and  the  ambient  steady  flow 
loss  terms  will  be  combined  to  form  a  simpler  expression 
below.  The  loss  mechanism  of  ambient  thermal  conduction 
through  the  gas  mixture  and  solid  stack  material  are  given  by 
^mixidTofdz)  and  K^oUdidToldz),  respectively,  while 
K]i^xud(dTQ/dz)  represents  an  additional  loss  by  thermal  con¬ 
duction  through  the  liquid  layer.  Use  of  Eqs.  (9)  and  (10) 
leads  to  considerable  simplifications  as  various  terms  com¬ 
bine. 
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The  acoustically  fluctuating  enthalpy  for  an  ideal  mix¬ 
ture  is  written  in  terms  of  temperature  and  concentration 
fluctuations,  and  is  given  by  w^=CpTQT*  +  »  where 

Cp  is  the  mixture  isobaric  specific  heat  and  is  the  mixture 
chemical  potential.  The  integration  over  the  solid,  liquid,  and 
gaseous  parts  of  the  stack  is  performed  and  the  expressions 
for  the  normalized  temperature,  axial  acoustic  velocity,  and 
concentration  fluctuations  are  substituted  and  time  averaged. 
With  these  modifications,  Eq.  (20)  becomes 

1  ^  [  (  Po  0  .  ^0  0 

+  Po^^l  Pc*”  ^^inix+ 


X 


Jo 


Irjdr^ 


Po  /  0  0  \ 

^  (^mix  ^liquid) 


PqPyI  ^2 


dTp 

^ liquid^  liquid*^  ^solid^Sf^dJ  ’ 


(21) 


where  the  integral  over  the  mixture  has  been  converted  to  an 
integral  over  a  single  pore  in  terms  of  the  dimensionless 
radial  coordinate  77.  Use  of  standard  thermodynamic  identi¬ 
ties  yields  further  simplifications 


H'c~^  M'^+  ^  ^  7’o(  •Suqdd)  >  (22a) 

Pi  Pi  Pi 


steady  flow  of  vapor  and  its  subsequent  return  flow  as  a 
liquid. 

Equation  (23)  is  the  most  convenient  form  for  evaluation 
from  first-order  acoustic  quantities,  but  is  not  the  best  form 
for  understanding  the  physical  properties  of  the  system.  If 
the  pqUW  term  in  Eq.  (20)  is  expressed  in  terms  of 
we  find 

PoC,  u*CpTpf*2vdv 

dcf^ 

■b  ro('^inbt~^liqiiid)^vapor“^inixro(^l~'y2)P0^12 

dTo 

(^mix^mix*^  ^liquid^  liquid ^soUd^  solid)  •  (24) 

This  form  resolves  the  additional  transport  terms  into  en¬ 
tropy  transport  due  to  the  mass  transport  of  the  mixture  and 
the  return  flow  of  the  liquid  and  entropy  transport  due  to  the 
ambient  diffusion  of  the  inert  gas  and  vapor  components. 

Evaluation  in  terms  of  the  volume  velocity  yields 


=  iRe[PAct^voi]+ 1  Re[^Act^voi(E,(kr) 


-beoEi(Xo))]- 


Cppo  dTQ 

t/vol 

2taA^  dz 

X  Im[E2(Xr) + eo<pE2(Xx))]  -  Gloss . 


(25) 


where 


and 


(M'Siix“’*'Squid)^~(^l~^2)~^^o(^2  ■^Uquid)- 

Pi  Pi 


(22b) 


where  xjj,,  x®.  ^Uquid  ^  the  ambient  entropies  per  unit 
mass  of  the  mixture,  vapor,  and  liquid  at  the  temperature  Tp , 
respectively.  With  the  identity  for  dC^ldz  from  Ref.  1,  Eq. 
(21)  becomes 

^AnyxRejpoC,  «*c^rof*277dj;j 
-  :5-AmixRe[^c,7'o(s2iix“'^Kquid)  f  u*Cll2ijdr) 

Z  Pi 
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■y  -o”o,  0  ,0 
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dz 
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[^ouxAinix"b /C|iqui(jA|iq,jj(j*l" *  (25) 


The  first  two  terms  represent  the  acoustically  transported  en¬ 
thalpy  Qkc+^ac-  The  first  term  can  be  evaluated  much 
like  Eq.  (11)  or  can  be  recognized  as  the  acoustic  heat  flux 
term  of  inert  gas  thermoacoustics.^-^  The  second  term  is 
To(5Siix“'^Uqind)  times  the  acoustic  mass  transport  of  vapor. 
The  integrals  for  acoustic  mass  transport  are  evaluated  in  Eq. 
(12).  The  third  term  is  the  enthalpy  loss  term  due  to  the 


_,o  X 

Cp  Po'^i 


(26) 


dTp 


Qloss  (AoiixAuijjf'l' A5Q|i(j/C5Qjjj“l*A|yq|,j^Kjjqujj) 


*1*  Aniii 
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<f>  P2^  (^  2  “  ■^liquid)^  12 


dTp 

dz 


,  (27) 


and  the  functions  Ei(Xy)  and  E2(X;)  are  defined  in  Eq.  (14) 
and  Eq.  (15). 

We  assume  that  the  liquid  layer  coating  all  surface  area 
of  the  stack  re{vesents  a  small  fraction  of  the  total  cross- 
sectional  area  of  the  stack  and  that  the  heat  conduction 
through  the  liquid  film  may  be  neglected.  This  assumption 
must  be  investigated  experimentally.  Since  the  axial  acoustic 
density  weighted  spedes  difference  velocity  is  shown  to  be 
zero  in  Ref.  1,  the  time-averaged  acoustic  power  is 


l^AC=5Re[PAct/voi]-  (28) 

This  is  the  first  term  in  Eq.  (25).  The  remainder  of  Eq.  (25)  is 
the  total  heat  flux  Q. 

Equations  (25)  and  (28)  provide  a  complete  description 
of  energy  fluxes  in  a  thermoacoustic  stack  with  an  inert  gas- 
condensing  vapor  working  fluid  in  a  pore  with  wet  walls 
under  the  approximations  described  above. 
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FIG.  1.  Mass  transport  directions  in  an  inert  gas-condensing  vapor  refrig¬ 
erator  for  different  values  of  the  temp^atiire  gradient  con^>med  to  the  heat 
transport  in  an  inert  gas  refrigerator  (2r  ^  thermal  hek  tr^port  flux, 
is  the  heat  transport  by  evaporation-condensation,  and  is  the 

mass  flux  of  vapor. 


V.  DISCUSSION 

Although  Eq.  (25)  is  quite  complex,  the  behavior  of  the 
different  terms  can  be  explained  in  terms  of  the  well- 
understood  inert  gas  thermoacoustic  refrigerator.  The  terms 
in  Eq.  (25)  containing  Xj-  arc  identical  to  those  in  Ref.  3  for 
inert  gas  refrigeration,  as  are  the  gas  and  solid  material  heat 
conduction  terms.  The  new  terms,  those  containing  X©  and 
the  additional  term  due  to  steady  flow  and  diffusion,  contain 
the  same  functions  evaluated  for  slightly  different  variables. 
Nfr  and  Nsc  are  both  on  the  order  of  1,  so  the  behavior  of  the 
heat  diffusion  and  mass  difiusion  terms  are  similar. 

Rgure  1  shows  that  in  an  inert  gas  refrigerator  with  the 
stack  properly  positioned,  the  acoustic  power  difference 
across  the  stack  is  converted  into  heat  flow  toward  the  hot 
end  of  the  stack  (maintained  near  room  temperature  by  a  heat 
exchanger).  This  removes  heat  from  the  cold  end  of  the  stack 
leading  to  cooling  below  room  temperature.  The  acoustic 
heat  flow  is  opposed  by  heat  conduction  in  the  gas  and  stack 
material  and  by  the  second  acoustically  driven  heat  term.  If 
losses  are  ignored,  cooling  is  possible  until  the  temperature 
gradient  is  large  enough  so  that  the  second  acoustic  heat  flux 
term  is  equal  to  and  opposite  in  direction  to  the  first  term. 
This  establishes  a  critical  gradient  for  refrigeration  and  that 
the  term  iRe[PAct^voi«i(^r)]  produces  a  heat  flux  away 
from  the  cold  end  of  the  stack  and  the  term  containing 
Im[H2(X7-)]  produces  a  heat  flux  toward  the  cold  end. 

Equation  (13)  for  the  mass  flux  contains  the  same  func¬ 
tions  as  the  heat  flux  but  evaluated  for  mass  diffusion  rather 
than  heat  diffusion.  The  vapor  flux  terms  are  in  the  same 


direction  as  the  heat  flux  terms.  When  the  temperature  gra¬ 
dient  is  small,  the  total  vapor  mass  flux  resulting  from  the 
two  diffusion  terms  is  in  the  same  direction  as  the  acoustic 
heat  transfer.  The  ideal  critical  gradient  for  the  vapor  transfer 
to  stop  and  reverse  directions  is  equal  to  the  critical  gradient 
for  heat  diffusion  multiplied  by  the  factor  I/9.  (p  is  usually 
greater  than  1.  Here,  “ideal”  means  that  the  third  term  in  Eq. 
(13)  describing  ambient  diffusion  can  be  neglected. 

The  vapor  transport  heat  flux  terms  in  Eq.  (25)  display 
the  same  behavior  as  the  mass  flux.  We  expect  ep  to  be 
positive,  since  the  entropy  of  the  mixture  should  be  greater 
than  the  liquid  entropy,  so  that  the  vapor  transport  terms  add 
to  the  acoustic  heat  transfer  terms  when  the  normalized  tem¬ 
perature  gradient  is  less  than  l/q>.  When  the  temperature  gra¬ 
dient  exceeds  the  critical  gradient  for  vapor  transport,  the 
flow  of  vapor  reverses  and  the  heat  transport  by  acoustic  heat 
difiusion  and  vapor  diffusion  oppose  each  other. 


VI.  IDEAL  COEFFICIENT  OF  PERFORMANCE  AND 
HEAT  TRANSFER 

A.  Theoretical  expressions 

The  coefficient  of  performance  (COP)  is  a  measure  of 
the  efficiency  of  a  refrigerator  and  is  equal  to  the  amount  of 
heat  pumped  from  a  cold  reservoir  divided  by  the  net  work 
input.  Thus,  the  coefficient  of  performance  is 


COP= 


AW 


(29) 


The  COP  can  be  evaluated  numoically  by  integrating 
through  the  stack  using  Eqs.  (1)  and  (2)  and  the  constancy  of 
Eq.  (25).  Then,  the  net  work  can  be  evalua^  from  Eq.  (28) 
evaluated  at  each  aid  of  the  stack  and  by  Qt^c  evaluated  at 
the  cold  end  of  the  stack  using  Eqs.  (25)  and  (28).  For  this 
preliminary  paper,  the  inert  gas-condensing  vapor  COP  is 
compared  to  the  inert  gas  COP  in  an  idealized  form  to  pro¬ 
vide  an  estimate  of  the  potential  improvements  provided  by 
the  wet-walled  system.  We  shall  consider  the  coefficient  of 
performance  of  an  inert  gas-vapor  thermoacoustic  refrigera¬ 
tor  in  the  inviscid,  boundary  layer,  and  standing  wave  limit. 

-r- 

The  amount  of  acoustic  woric  used  to  pump  Qc  is  the  differ¬ 
ence  in  acoustic  intensity  on  the  hot  and  cold  sides  of  the 
stack.  Thus, 


■Amix[(^AC^)right  (f^AC^)left] 


d 


=  -“|^(Re[PAC^voi]) 


=  ->:rRe^P 


AC 


dUyoi  dP AC  - 
+  - ; - 


dz 


dz 


(30) 


where  we  have  used  the  fact  that  ^ 

stack  length.  The  coordinates  are  indicated  in  Fig.  1.  This 
approximation  puts  an  upper  bound  on  COP  since  it  assumes 
losses  in  other  parts  of  the  refrigerator  are  negligible. 
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In  the  standing  wave  limit,  f/voi=£^voi,s  and  Pac 
=  iTac.s  •  Th®  boundary  layer  limit  of  the  f  (Xy)  functions  is 


F(,kj)=  l-il-i)Sj/R,  (31) 


with  j  taken  as  T  referring  to  thermal  effects  or  D  referring  to 
mass  difiusion  effects. 

Applying  these  limits  to  the  change  in  acoustic  work 
yields 


£iW=- 


L 

^2 

\St  y 

dp 

2 

[R  n^l  y-1 

R 

P AC,s^voi,5  dTp  [  Sr  ^  ^2  y  Sp 

To  dz  [R  ^ y—l  R 


(32) 


The  acoustic  heat  pumped  from  the  cold  reservoir  can  be 
written  with  these  limits  as 


;r 

Qc  2 


^£>1 


^ppo  r,2 


Sf 

-j^+eo'P 


(33) 


Taking  the  ratio  of  Eq.  (33)  and  Eq.  (32)  and  factoring  out 
”“^ac;5^voi,5/2  from  the  numerator  and  —L(oA^ 
X  from  the  denominator  yields 


TOP- 


{Sr+eM- 

- — - — "7 - H - 


(34) 


Two  factors  used  in  the  literature^  of  inert  gas  thermoacous¬ 
tics  can  be  identified: 

the  dry  critical  temperature  gradient 


^u>A^Pac;s 

c^Pof/vou’ 

and  the  normalized  temperature  gradient 


(35) 


"^T^^dTpIdz 

^^crit  ^^crit 


(36) 


The  critical  temperature  gradient  is  the  ambient  temperature 
gradient  at  which  the  temperature  oscillations  in  the  mixture 
due  to  acoustic  pressure  changes  are  the  same  as  the  local 
stack  temperature  when  the  acoustic  wave  displaces  the  mix¬ 
ture  along  the  stack.  In  other  words,  at  the  critical  gradient  a 
parcel  of  gas  in  the  mixture  will  always  be  in  thermal  equi¬ 
librium  with  the  local  stack  temperature  and  thus  will  not 
exchange  thermal  energy  with  the  stack.  However,  in  a  wet 
system  heat  is  transported  by  the  vapor,  as  long  as  the 
changes  in  the  partial  pressure  of  vapor  due  to  acoustic  pres¬ 
sure  changes  are  different  from  the  vapor  pressure  at  the 
stack  wall.  The  partial  pressure  at  the  wall  depends  only  on 
the  temperature  and  the  latent  heat  through  the  parameter  <p. 
We  further  identify  and  approximate  the  Carnot  coeffi¬ 
cient  of  performance  as 


COPcamot 


Tc  _  To  dz 

Th-Tc^  LdTo^ 


(37) 


This  approximation  assumes  that  the  temperature  gradient  is 
linear  and  that  the  cold  temperature  and  the  stack  average 
temperature  are  comparable.  Equation  (34)  can  then  be  ex¬ 
pressed  as 
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(38) 


If  there  is  no  vapor  in  the  system,  n2=sn==0  and  the  coef¬ 
ficient  of  performance  reduces  to  the  expression  COP 
=  rCOPcanjot,  which  agrees  with  Swift’s  expression  for  inert 
gas  thermoacoustics.^ 

This  idealized  expression  is  a  valuable  aid  to  under¬ 
standing  the  transport  properties.  The  relationship  of  the 
change  in  sign  of  the  mass  difiusion  acoustic  transport  as  F 
exceeds  l/(p  is  clear.  The  ideal  COP  will  be  increased  for 
systems  in  which  the  normalized  temperature  gradient  is  be¬ 
tween  zero  and  I/9.  If  the  normalized  gradient  exceeds  l/<p, 
the  mass  difiusion  effects  oppose  the  heat  diffusion  effects.  If 
«Z)^(7^y~l)('*2^'*i)  Ihe  normalized  gradient  is  be¬ 
tween  0  and  l/<p,  the  COP  of  the  wet  system  will  be  greater 
than  the  COP  of  the  dry  system  with  the  same  dry  properties. 


B.  Example  analysis  for  selected  inert  gas-vapor 
mixtures 

As  an  indicator  of  the  potential  of  inert  gas-condensing 
vapor  mixtures  to  improve  thermoacoustic  refrigeration,  we 
have  calculated  the  ideal  COPs  relative  to  Carnot  from  Eq. 
(38)  for  a  thermoacoustic  cooler  operating  between  280  and 
310  K  as  a  function  of  the  normalized  temperature  gradient 
r.  Physical  properties  of  the  pure  gas  and  vapor  at  standard 
state  and  as  functions  of  temperature  may  be  found  in  nu¬ 
merous  reference  books. Transport  properties,  such  as 
viscosity,  thermal  conductivity,  or  diffusion  coefficient,  of 
pure  gases  and  binary  mixtures  of  gases  are  calculated  based 
on  gas  dynamic  theory.*^  Thermophysical  properties  such  as 
the  specific  heat  and  the  ratio  of  specific  heats  of  the  mixture 
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FIG.  2.  Coefficient  of  performance  relative  to  the  Carnot  coefficient  of 
performance  versus  normalized  temperature  gradient  for  pure  helium 

_ helium  and  water  vapor - ,  helium  and  ammonia  *  *  *  *,  helium 

and  butane  -*-♦  helium  and  R134a  — . 


FIG.  3.  Normalized  cooling  power  versus  normalized  temperature  gradient 

for  pure  helium - helium  and  water  vapor - ,  helium  and  ammonia 

•  •  *  helium  and  butane  — ,  helium  and  R134a  — . 


are  calculated  with  mole  firaction  weighting.^^  The  mole  fiac- 
tions  of  gas  and  vapor  in  the  mixture  are  calculated  from  the 
Clausius-Clapeiyon  equation  as  a  function  of  ambient  local 
temperature.^"^  When  calculating  the  mole  fractions  it  is  as¬ 
sumed  that  the  mean  pressure  in  the  model  refrigerator  is  a 
constant;  thus,  an  air  and  water  vapor  mixture  would  have 
equal  mole  fractions  for  a  mean  pressure  of  2  atmospheres  at 
373  K.  The  numerator  of  the  coefficient  of  performance  rela¬ 
tive  to  Camot,  Eq.  (38),  is  the  heat  pumped  from  the  cold 
side  of  the  stack.  Thus,  the  thermophysical  quantities  of  the 
mixture  in  the  numerator  are  calculated  at  an  ambient  tem¬ 
perature  of  280  K.  The  denominator  of  Eq.  (38)  is  the  work 
required  to  pump  the  heat  and  is  evaluated  at  the  average 
temperature  across  the  stack  (295  K). 

Figure  2  indicates  that  the  coefficient  of  performance  of 
the  helium-condensing  vapor  mixture  is  higher  than  that  of 
helium  for  low  normalized  temperature  gradients.  The  pure 
helium  system  can  achieve  higher  COPs  relative  to  Camot  at 
higher  critical  gradients,  but  as  will  be  demonstrated  below, 
for  much  lower  cooling  powers  than  the  helium-ammonia, 
helium-butane,  or  helium-R134a  mixtures. 

Swift^  suggests  that  the  potential  for  different  gas  mix¬ 
tures  to  transport  energy  can  be  investigated  by  assuming 
|PacHp^|wacI  that  the  highest  acoustic  Mach  number 
(M=|wac/^I)  achievable  in  an  acoustic  system  is  constant 
at  approximately  0.1.  Further,  we  set  the  ratio  SjlR  for  each 
system  equal.  This  corresponds  to  optimizing  the  stack  spac¬ 
ing  for  each  mixture.  Under  these  assumptions  Eq,  (33)  be¬ 
comes 

lecl=^^(y)[(i-n+Mi-«^r)].  (39) 

In  Fig.  3  we  plot  the  normalized  cooling  power 

ipc\ciium 


I  (2  i  normalized 


The  relative  cooling  power  for  the  pure  helium  system  is 
then  given  by  (l-F).  The  thermophysical  properties  for  this 
calculation  are  evaluated  at  10  atmospheres  and  280  K. 
Large  gains  in  heat  flux  from  the  cold  heat  exchanger  are 
predicted  at  relatively  high  COPRs.  These  results  are  encour¬ 
aging,  but  a  more  realistic  calculation  should  include  the 
additional  stack  length  and  attendant  increased  losses.  Initial 
estimates  of  the  new  loss  terms  in  Eq.  (25)  which  have  not 
been  included  in  this  analysis  indicate  that  they  are  at  most 
on  the  order  of  the  heat  conduction  losses  in  inert  gas  ther¬ 
moacoustics,  but  this  must  also  be  confirmed  by  numerical 
calculations. 

VIL  CONCLUSIONS 

The  transport  equations  for  wet  thermoacoustics  have 
been  developed.  The  heat  energy  transported  and  the  coeffi¬ 
cient  of  performance  are  both  increased  for  wet  stacks  oper¬ 
ated  between  a  temperature  gradient  of  zero  and  the  critical 
temperature  gradient  for  mass  transport.  The  cooling  by  va¬ 
por  transport  is  analogous  to  a  conventional  vapor  cycle  re¬ 
frigerator.  High  entropy  vapor  is  transported  from  the  cold 
side  of  the  stack  to  the  ambient  where  it  condenses.  This 
system  combines  thermoacoustic  refrigeration  with  vapor 
cycle  refrigeration  in  a  simple  acoustically  driven  device. 

Several  experimental  and  calculational  questions  must 
be  answered  before  it  can  be  determined  if  inert  gas- 
condensing  vapor  thermoacoustics  can  be  exploited  to  im¬ 
prove  thermoacoustic  refrigerators. 

(1)  A  full  numerical  calculation  must  be  performed  to  see  if 
the  improvement  in  cooling  power  shown  above  can  be 
achieved  in  a  realistic  refrigerator.  In  this  paper  we  have 
derived  the  necessary  relations  for  calculations  in 
DeltaE.* 

(2)  The  behavior  of  the  condensed  liquid  layer  in  the  stacks 
must  be  investigated  experimentally  for  a  variety  of 
stack  materials  and  vapor  substances.  We  have  per¬ 
formed  our  calculations  with  the  simplest  assumption: 
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that  the  liquid  layer  is  thin  everywhere  and  that  natural 
flow  occurs  in  the  steady  state  to  return  the  transported 
and  condensed  vapor. 

(3)  The  behavior  of  the  heat  and  mass  transfer  at  the  heat 
exchangers  must  be  investigated.  Our  analysis  of  en¬ 
thalpy  flow  is  only  valid  within  the  stacks. 

(4)  The  assumption  that  the  temperature-driven  diffusion  ef¬ 
fect  (Jtj-  terms)  are  negligible  should  be  investigated  in 
more  detail.  We  have  relied  on  the  analysis  of  Ref.  6  in 
our  neglect  of  jfcr ,  but  kj  may  be  larger  in  a  given  can¬ 
didate  system. 

There  are,  of  course,  other  questions  and  considerations  be¬ 
sides  those  listed  above.  Nonetheless,  this  is  a  promising  and 
interesting  field  of  investigation. 
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(57)  ABSTRACT 

A  thermoacouslic  refrigeration  device  employs  a  gas-vapor 
mixture  as  the  working  fluid.  As  a  result,  the  refrigeration 
device  operates  according  to  a  modified  thermoacoustic 
refrigeration  cycle  that  adds  a  condensation-vaporization 
cycle  to  the  thermoacoustic  cycle.  The  resulting  modified 
refrigeration  cycle  increases  the  efficiency  of  heat  transport 
by  harnessing  the  translational  motion  of  the  vapor,  as  well 
as  the  usual  acoustic  oscillations,  to  transport  the  heat  energy 
from  one  end  of  a  thermal  stack  to  the  other. 

10  Claims,  2  Drawing  Sheets 
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THERMOACOUSTIC  REFRIGERATION 
DEVICE  AND  METHOD 

This  application  claims  the  benefit  of  provisional  appli¬ 
cation  No.  60/334,945,  filed  Dec.  4,  2001.  5 

BACKGROUND  OF  THE  INVENTION 

1.  Field  of  the  Invention 

This  invention  relates  to  improvements  in  thermoacoustic 
refrigeration,  and  in  particular  to  a  thermoacoustic  refrig-  IC 
eration  device  which  employs  a  gas-vapor  mixture  as  the 
working  fluid. 

The  invention  also  relates  to  a  refrigeration  method,  and 
in  particular  to  a  thermoacoustic  refrigeration  method  that 
modifies  the  conventional  ‘thermoacoustic  refrigeration  15 
cycle  by  adding  condensation  and  vaporization  to  the  ther¬ 
moacoustic  cycle. 

In  conventional  thermoacoustic  refrigeration,  heat  energy 
is  transported  primarily  or  solely  by  waves  acoustically 
induced  in  an  inert  gas.  The  waves  may  be  standing  waves  20 
or  traveling  waves  that  oscillate  or  travel  from  one  side  of 
the  thermal  stack  to  the  other,  but  in  either  case  are  arranged 
to  exploit  the  temperature  differential  between  areas  of 
compression  and  areas  of  rarefaction  relative  to  a  thermal 
stack  or  regenerator,  the  thermal  stack  transferring  heat  25 
energy  from  the  cold  areas  to  the  hot  areas  to  achieve 
refrigeration.  The  purpose  of  the  vaporization-condensation 
cycle  of  the  invention  is  to  increase  the  efficiency  of  heat 
transport  by  harnessing  the  translational  motion  of  the  vapor, 
as  well  as  the  usual  acoustic  oscillations,  to  transport  the 
heat  energy  from  one  end  of  the  thermal  stack  to  the  other, 
i.e.,  by  using  acoustic  mass  transfer  as  well  as  acoustic  heat 
transfer  to  transport  heat  energy  up  the  heat-absorbing  stack. 

2.  Description  of  Related  Art 

The  basic  principles  of  thermoacoustic  refrigeration  have 
been  known  for  more  than  a  decade.  In  its  most  basic  form, 
thermoacoustic  refrigeration  is  a  process  that  utilizes  acous¬ 
tic  energy  to  pump  or  transport  heat  through  a  thermal  stack 
between  a  cold  heat  exchanger  at  one  end  and  a  hot  heat 
exchanger  at  the  other  end,  the  acoustic  energy  being  in  the 
form  of  standing  or  traveling  acoustic  waves  generated  by  a 
loudspeaker  or  similar  moving  part  to  cause  the  mechanical 
compression  and  expansion  of  a  working  fluid  needed  for 
the  cooling  cycle. 

In  all  such  devices,  heat  is  transported  solely  by  acoustic 
waves  in  an  inert  gas.  Since  the  acoustic  wave  generator  is  45 
the  only  moving  part  in  a  thermoacoustic  refrigeration 
device,  thermoacoustic  refrigeration  devices  have  the  poten¬ 
tial  for  greater  reliability,  smaller  size,  and  lower  cost  than 
conventional  refrigeration  devices.  Despite  these  potential 
advantages,  however,  thermoacoustic  refrigeration  has  one  50 
substantial  disadvantage  that  has  heretofore  prevented  more 
widespread  application  of  the  technology,  namely  relatively 
low  cooling  power  in  comparison  with  conventional  vapor 
cycle  refrigeration  devices.  Essentially,  the  problem  is  that 
the  oscillations  in  temperature  resulting  from  application  of 
an  acoustic  wave  to  the  inert  gas  used  as  a  worldng  fluid  are 
relatively  small  and  thus  have  a  limited  ability  to  transfer 
heat. 

The  basic  standing  wave  version  of  the  conventional 
thermoacoustic  device,  illustrated  in  FIG.  1,  is  simply  a 
hollow  resonance  tube  1  filled  with  an  ordinary  inert  gas  and 
having  a  speaker  or  acoustic  driver  2  at  one  end,  a  hard 
termination  3  at  the  other  end,  and  a  dry  stack  or  regenerator 
4  sandwiched  between  hot  and  cold  heat  exchangers  5,6.  The 
dry  stack  or  regenerator  is  composed  of  a  solid  material  and 
finely  divided  into  sections  or  passages  with  which  the  65 
working  fluid  exchanges  heat  with  the  solid  material.  In  the 
case  of  a  standing  wave,  a  temperature  differeniial  is  estab- 


2 

lished  between  areas  of  compression,  in  which  heat  is 
transferred  from  the  working  fluid  to  the  thermal  stack,  and 
rarefaction,  in  which  heat  is  transferred  from  the  thermal 
stack  to  the  working  fluid. 

Details  of  the  refrigeration' device  shown  in  FIG.  1  may  be 
found  in  T.  Hofler,  "^Thermoacoustic  refrigeration  design 
and  performance^  Ph.D.  dissertation.  Physics  Department, 
University  of  California  at  San  Diego,  1986,  while  addi¬ 
tional  information  and  background  on  standing  or  traveling 
wave  refrigeration  devices  may  be  found  in  U.S.  Pat,  Nos. 
4,114,380  (Ceperley);  43,555,517  (Ceperley);  4,398,398 
(Wheatley  et  al.);  4,489,553  (Wheatley  et  al.);  4,722,201 
(Hofler  et  al.);  4,858,441  (Wheatley  el  al.);  4.953,366  (Swift 
et  al.);  5,165,243  (Bennett);  5,647,216  (Garrett);  5,673,561 
(Moss);  5,901,556  (Hofler);  6,032.464  (Swift  et  al.);  6,164, 
073  (Swift  et  al.);  and  6,233,946  (Matsuda). 

The  present  invention,  in  contrast,  augments  the  purely 
acoustic  heal  transfer  of  previously  proposed  thermoacous- 
lic  refrigeration  devices  by  adding  an  evaporation- 
condensation  cycle  similar  to  that  of  conventional  vapor- 
cycle  based  refrigeration  devices,  but  without  the 
mechanical  complexity  of  the  conventional  refrigeration 
device.  This  is  achieved  by  mixing  vapor  with  an  appropri¬ 
ate  gas  working  fluid,  and  by  permitting  the  working  fluid  to 
locally  evaporate  and  condense  on  the  hot  and  cold  ends  or 
sides  of  the  thermal  stack,  thereby  transferring  heat  energy 
by  acoustic  mass  transfer  as  well  as  by  acoustic  heat  transfer. 

With  the  proper  selection  of  the  gas-vapor  mixtures  used 
as  working  fluids  in  the  thermoacoustic  re^geration  devices 
of  the  invention,  criteria  for  which  are  described  below, 
predicted  heat  pumping  power  and  coefficient  of  perfor¬ 
mance  relative  to  the  Carnot  cycle  can  be  increased  for  an 
inert  gas-vapor  working  fluid  compared  to  a  similar  purely 
inert  gas  woildng  fluid. 

SUMMARY  OF  THE  INVENTION 

It  is  accordingly  an  objective  of  the  invention  to  provide 
a  thermoacoustic  refrigeration  device  that  provides 
improved  performance  without  a  substantial  increase  in 
complexity. 

It  is  also  an  objective  of  the  invention  to  provide  a 
thermoacoustic  refrigeration  method  that  increases  the  heat 
transfer  capability  of  the  thermoacoustic  refrigeration  device 
by  modifying  the  conventional  thermoacoustic  refrigeration 
cycle  to  include  evaporation  and  condensation,  thereby 
enlisting  translation  of  vapors  as  a  mechanism  for  heat 
transfer  in  addition  to  the  conventional  acoustic  heat  trans¬ 
fer. 

It  is  a  still  further  objective  of  the  invention  to  provide 
criteria  for  selecting  a  gas-vapor  mixture  that  can  be  used  in 
a  thermoacoustic  refrigeration  .device  in  order  to  increase 
performance  by  making  use  of  translating  vapor  to  cany 
heal  energy,  and  condensation  of  the  vapor  to  transfer  the 
heat  energy  to  the  refrigeration  stack. 

According  to  a  preferred  embodiment  of  the  invention, 
the  thermoacoustic  refrigeration  device  of  the  invention  may 
include,  as  in  a  conventional  thermoacoustic  refrigeration 
device,  a  resonance  tube  filled  with  a  working  fluid,  an 
acoustic  driver  at  one  end,  a  hard  termination  at  the  other 
end,  a  thermal  stack  situated  within  the  tube  and  made  up  of 
a  finely  divided  structure  composed  of  a  solid  material,  and 
cold  and  hot  heat  exchangers  at  opposite  ends  of  the  stack. 
Unlike  the  inert  gas  of  the  conventional  device,  however,  the 
working  fluid  employed  by  the  refrigeration  device  of  the 
preferred  embodiment  is  a  gas-vapor  mixture,  the  stack  is 
composed  of  a  solid  material  that  is  wettable  by  condensed 
vapor,  and  means  are  provided  to  return  condensed  vapor 
from  the  cold  side  of  the  stack  to  the  hoi  side.  The  return 
means  may  be  a  pump  or,  for  simplicity,  a  wick. 
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Alternatively,  the  stack  may  be  arranged  to  permit  return 
solely  or  primarily  by  gravity. 

The  method  of  the  preferred  embodiment  of  the  invention 
modifies  conventionsJ  thermoacoustic  cooling,  which  sim¬ 
ply  involves  generating  acoustic  waves  in  the  working  fluid, 
by  modifying  the  thermoacoustic  refrigeration  cycle  to 
include  at  least  the  following  steps: 

a.  In  response  to  an  applied  acoustic  wave,  a  parcel  of  gas 
in  the  working  fluid  is  caused  to  undergo  translation 
along  the  stack  and  consequent  acoustic  compression, 
thereby  decreasing  the  parcel’s  volume  and  increasing 
its  temperature; 

b.  The  decreased  volume  and  increased  temperature 
increases  the  partial  pressure  of  the  vapor  within  the 
parcel; 

c.  The  parcel  then  slows,  stops,  and  reverses  its  transla¬ 
tional  motion,  while  at  the  same  time  exchanging  heat 
and  vapor  with  the  stack  as  a  result  of  the  parcel’s 
increased  temperature  relative  to  the  stack; 

d.  At  the  time  of  reversal,  the  increased  partial  pressure 
relative  to  the  vapor  pressure  at  the  stack  wall  causes 
vapor  to  condense  from  the  parcel  to  the  stack  plate; 

e.  The  gas  parcel  then  undergoes  acoustic  rarefaction  and 
is  translated  back  past  the  ambient  position,  increasing 
its  volume  and  decreasing  its  temperature; 

f.  The  acoustic  rarefaction  in  tum  causes  a  decrease  in 
partial  pressure  of  vapor  within  the  parcel; 

g.  The  parcel  again  slows,  stops  and  reverses  its  transla¬ 
tional  motion  while  exchanging  heat  and  vapor  with  the 
stack,  this  time  absorbing  heat  from  the  stack;  * 

h.  Since  the  partial  pressure  of  the  vapor  in  the  parcel  is 
lower  than  the  vapor  pressure  at  the  stack  wall,  the 
vapor  will  evaporate  to  the  parcel  from  the  liquid  layer 
coating  the  stack. 

Because  of  the  modifications  to  the  conventional  iher- 
moacoustic  refrigeration  cycle,  heat  is  transported  from  one 
end  of  the  stack  to  the  other,  both  as  a  result  of  the 
temperature  differentials  resulting  from  compression  and 
rarefaction  of  the  gas,  and  as  a  result  of  vapor  moving  up  the 
stack  and  exchanging  heat  at  each  end. 

The  performance  of  the  refrigeration  device  of  the  pre¬ 
ferred  embodiment  depends  on  the  characteristics  of  the 
gas-vapor  mixture.  According  to  another  aspect  of  the 
invention,  the  gas-vapor  mixture  is  chosen  to: 


per  mole,  T^  is  the  ambient  temperature,  and  is  the 
universal  gas  constant.  This  parameter  controls  the  critical 
gradient  and  therefore  the  stack  length  necessai^  to  produce 
a  given  temperature  difference. 

5  Although  a  particular  standing  wave  device  is  illustrated 
herein,  those  skilled  in  the  art  will  appreciate  that  the  device 
and  the  method  of  the  invention  may  also  be  applied  to 
traveling  wave  refrigeration  techniques,  and  that  the  simple 
tube,  stack,  and  heat  exchange  structure  of  the  illustrated 

10  embodiment  may  be  freely  varied  by  those  skilled  in  the  art 
in  accordance  with  the  principles  described  in  any  of  the 
above  references,  without  departing  from  the  scope  of  the 
invention,  so  long  as  the  modified  device  employs  a  gas- 
vapor  mixture  as  the  working  fluid. 


15  BRIEF  DESCRIPTION  OF  THE  DRAWINGS 

RG.  la  is  a  schematic  diagram  of  a  conventional  standing 
wave-type  thermoacoustic  refrigeration  device. 

RG.  lb  is  a  schematic  diagram  of  a  thermoacoustic 
refrigeration  device  constructed  in  accordance  with  the 
principles  of  a  preferred  embodiment  of  the  invention. 

RG.  2a  is  a  graph  of  the  ideal  coefficient  of  performance 
relative  to  Carnot  plotted  against  the  normalized  temperature 
gradient, 

RG.  2b  is  a  graph  of  the  normalized  cooling  power  for  the 
inert  gas- vapor  mixtures  whose  coefficient  of  perfoimance  is 
plotted  in  RG.  2a, 

DETAILED  DESCRIPTION  OF  THE 
PREFERRED  EMBODIMENTS 

As  illustrated  in  RG.  lb,  which  is  a  conceptual  drawing 
of  a  refrigeration  device  constructed  according  to  the  prin¬ 
ciples  of  a  preferred  embodiment  of  the  invention,  the 
thermoacoustic  refrigeration  device  of  the  preferred  embodi¬ 
ment  may  include  a  resonance  tube  10  filled  with  a  working 
fluid,  an  acoustic’driver  11  at  one  end,  a  hard  termination  12 
at  the  other  end,  a  thermal  stack  13  situated  within  the  tube 
and  made  up  of  a  finely  divided  structure  composed  of  a 
solid  material,  and  respective  cold  and  hot  heat  exchangers 
14,15  at  opposite  ends  of  the  stack.  The  heat  exchangers  may 
include  conduits  not  shown)  through  which  are  supplied  hot 
and  cold  liquids,  or  any  other  heat  conducting  structure,  for 
the  purpose  of  transporting  heat  away  from  the  thermal 
stack,  in  the  manner  of  a  conventional  heat  exchanger.  The 
working  fluid  employed  by  the  refrigeration  device  of  the 
a.)  maximize  the  amount  of  heal  carried  by  mass  relative  45  preferred  embodiment  is  a  gas-vapor  mixture  whose  com- 
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to  the  amount  of  heat  carried  by  thermal  effects,  rep¬ 
resented  by  €£,  defined  by  the  relationship: 


1 

£0— - 

Cj,  Po  n\  y 


50 


where  is  the  mass  density  of  the  vapor  in  the  mixture, 
is  the  mass  density  of  the  mixture,  n^  is  the  number  density 
of  the  mixnire,  nj  ^  is  the  number  density  of  the  inert  gas,  s„,^ 
is  the  entropy  per  unit  mass  of  the  mixture  the  entropy 

per  unit  mass  of  the  condensed  liquid,  and  is  the  heat 

capacity  at  constant  pressure  per  unit  mass  of  the  mixture 
and, 

b.)  minimize  the  heat  transfer  coefficient  W  defined  by  the 
relationship: 

y  RoTo' 
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60 


65 


where  y  is  the  ratio  of  the  specific  heats  of  the  gas  and  of  the 
mixture;  1  is  the  latent  heat  of  vaporization  of  the  mixture 


position  meets  the  criteria  described  below,  while  the  stack 
is  composed  of  a  solid  material  that  is  wettable  by  condensed 
vapor  such  that  a  thin  sheet  of  liquid  will  cover  the  stack. 
Finally,  means  are  provided  to  rfetum  condensed  vapor  from 
the  cold  side  of  the  stack  to  the  hot  side.  The  return  means 
may  be  a  pump  16,  as  illustrated  or,  for  simplicity,  a  wick. 
Alternatively,  the  stack  may  be  arranged  to  permit  rerum 
solely  or  primarily  by  gravity. 

According  to  the  preferred  embodiment  of  the  invention, 
a  modified  thermoacoustic  refrigeration  cycle,  including 
vaporization  and  condensation  steps,  is  provided  as  follows: 

a.  In  response  to  an  acoustic  wave  applied  by  acoustic 
driver  11,  a  parcel  of  gas  in  the  working  fluid  is  caused 
to  undergo  translation  along  the  stack  and  consequent 
acoustic  compression,  thereby  decreasing  the  parcel’s 
volume  and  increasing  its  temperature; 

b.  The  decreased  volume  and  increased  temperature 
increases  the  partial  pressure  of  the  vapor  within  th^ 
parcel; 

c.  The  parcel  then  slows,  stops,  and  reverses  its  transla¬ 
tional  motion,  while  at  the  same  time  exchanging  heat 

vapor  with  the  stack  13  as  a  result  of  the" parcel’s 
increased  temperature  relative  to  the  stack; 
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d.  At  the  lime  of  reversal,  the  increased  partial  pressure 
relative  to  the  vapor  pressure  at  the  stack  waif  cau^s 
vapor  to  condense  from  the  parcel  to  the  adjacent  stack 

undergoes  acoustic  rarefaction  and 

ts  v^nmfinfa  P""'  increasing 

volume  and  decreasing  its  temperature; 

rarefaction  in  turn  causes  a  decrease  in 
^lal  pressure  of  vapor  within  the  parcel; 
o-  e  parcel  again  slows,  stops  and  reverses  its  transla 
exchanging  heat  and  vapor  wift  ie 
stack,  this  time  absorbing  heat  from  the  stack- 
h.  Since  Ae  partial  pressure  of  the  vapor  in  the  parcel  is 
lower  thM  the  vapor  pressure  at  the  stack  wall  the 

co^no''|S/s''S^  iWer 

theimodyn^c  cycle  describes  how  an  inert  »as- 

and  vapor  from  om 

stark  “  uhosfn  so  that  the 

stack  remains  wet  over  the  temperature  range  of  interest 
WiA  the  properly  chosen  inert  gas  and  con^nsino  vapor 
heat  ener^  is  earned  up  the  slack  by  acoustic  mLsIransfer  20 
as  well  as  by  acoustic  heat  transfer.  TOe  vapor  mo  Wng  Se 
stack  condenses  on  the  walls  of  the  slack  and  wSl 

Md  either  by  gravity  or  by  an  external  wicking  anSig? 
ment  or  pump  16.  ®  "'••uigc 

The  equations  describing  acoustic  propagation  in  an  inert  •>«  n,;..  •  .c 

gas-vapor  mixture  in  a  thermoacoustic  stack  are  the  acoustic  relative  coefficient  of  performance 

wave  equauon  and  the  heat  pumping  equation,  both  of  which  of  the  ^  5  different  working  fluids  as  a  function 

wLTSitfiT’  “ f described  here.  From  these  2^ temperature  gradient  is  depicted  in  FIG. 
equations,  figures  of  ment  such  as  the  coefficient  of  perfor-  a  '  ■  properties  of  the  working  fluids  such  as 

m^ce  and  the  normalized  cooling  power  may  be  dSv^"  "“"jber  density,  and  specific  heat,  were  c’alSat^ 

Si  ffie  moT  handbooks, 

^  r  pressure  is  assumed  to  be  a  constant  at  10 

effi^enc^rfa  fffrhf  perfonnMce  (COP)  is  a  measure  of  the  for  all  working  fluids.  Since  the  numerator  of  the 

pumped  from  a  amount  of  beat  't^ealtzed  efficiency  represents  the  heat  pumped  from  a  cold 

pumped  from  a  cold  reservoir  divided  by  the  net  work  input  reservoir  the  working  fluid  properties  in  the 

upper  bound  for  this  efficiency  mav  be  35  evaluated  at  a  temnpratiirA  in.-r  oor\  ip  tna 


where  p,"  is  the  mass  density  of  the  vapor  in  the  mixture  p 
he  m^s  densuy  of  the  mixture,  n„  is  the  number  density 
of  the  mixnire,  n,  is  the  number  density  of  the  inen  oas  s  ^ 
IS  the  entropy  per  unit  mass  of  the  mixie  S,^„„  the  en&oov 
10  per  umt  mass  of  the  condensed  liquid,  and  c  is  the  hei 
ca^city  at  constant  pressure  per  unit  mass  of  the  mixture 
ihe  expression  for  the  idealized  COPR  may  now  he 
considered  m  more  detail.  If  there  is  no  vapor  withinTaJ 
^ture  the  Co  and  n/  terms  go  to  zero  and  the  expipssfon 
reduces  to  the  accepted  idealized  efficiency  for  a  ffiy  ther- 

^  gas-vapor  mixture  that  the 
efficiency  will  be  mcreased  relative  to  a  comparable  drv 
imxture  for  values  of  the  normalized  temperature  oradi^ 
between  zero  and  I/4,  and  if  is  greater  tiian  ^ 


y  »! 
r-lnf 


(7°  C),  while  the 


^'Xwr  tound"!^  properties  i 

makinTLS  knV  ®®“ency  may  be  calculated  by  evaluated  at  a  temperature  of  280  K 

acousti^sSnrw^nfirasem-^  ^  ??22^cT  teiSp^e  of^^S 

y  made  by  assuming  a  linear  temnftrahir#»  nrn 


TTiav  a  runner  simplification 

•♦k'  assuming  a  linear  temperature  «T^dient 

efficiency 

1? thM  ™re  with  the  above  simplifications 


fr,  ^  *at  the  idealized  efficiency  is 

'forking  fluids  cLpared  to  a 
p  e  helium  working  fluid  at  low  values  of  the  normalized 
temperature  ^adient.  For  the  helium-butane  andTelium 
working  fluids  at  low  values  of  F  the  idealized 
efficiency  may  be  increased  by  a  factor  of  3.  FIG  2a  also 
demonstrates  the  dependence  of  the  COPR  on  the  value  of 
^  for  the  working  fluids  under  consideration.  The  COPR  of 
.  „ .  rnp^  gas-vapor  working  fluid  is  greater  than  the  dry 

where  F  is  the  actual  temperature  gradient  across  ihp  ei^pk  50  ^  as  the  normalized  temperature  ®radient  lies 


copR=r‘ 


~  F) + £o^£i(1 

y 


-ipV) 


nn^,r  A  u  t  T  giaoient  across  the  stack 

normalized  by  the  dry  cmical  temperature  gradient,  6.  and 

bo  are  the  thermal  and  difiiision  penetration  depths  within 
ffie  mixture  yis  the  ratio  of  specific  heats  for  tte  mixture. 
Hj  and  n,  are  the  number  densities  of  the  inert  vas  and 


55 


_  y-1  / 

^  7  KT,' 

Where  F  is  the  ratio  of  ±e  specific  heats  of  the  gas  and  of  the 
^xture;  1  is  t^he  latent  heat  of  vaporization  of  the  mixture  per 
mole,  T„  ,s  the  ambient  temperature,  and  R„  is  the  umvprsal 
<=as  constant.  An  expression  that  describes  the  relative 
^ount  of  heat  transported  by  thermal  or  mass  diffusion 


oustic  refngerator  is  an  expression  for  the  normalized 
cooling  power  given  by,  ‘iunnaiizea 

Helium 


60  vaZhtll  fluid  and  other 

anabJes  are  the  same  as  defined  before.  Standard  assump- 

lions  were  used  to  derive  this  formula  such  as  assuming  a 
Ihe  acousuc  Mach  number  and  assuming  the  ratio  of 
e  thermal^  penetration  depth  in  the  working  fluid  to  the 

65  Within  the  Stack  is  the  same 

for  the  wet  and  diy  working  fluids. 

The  noim^iz^  cooling  power  for  various  working  fluids 
IS  depicted  in  HG.  2i  as  a  ffinction  of  the  non^Sd 
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temperature  gradient.  Xhe  working  fluids  considered  in  this 
calculation  are  the  same  as  those  shown  in  FIG.  2a.  'Die 
ihennal  properties  of  the  working  fluids  are  calculated  using 
standard  formula  from  gas  property  handbooks  at  a  mean 
pressure  of  10  bar  and  an  ambient  temperature  of  280K  (7° 
C-).  The  normalized  cooling  power  for  the  helium-butane 
and  helium-R134a  working  fluids  at  low  values  of  F  is 
approximately  a  factor  of  3  greater  than  the  pure  helium 
working  fluid.  FIGS.  2a  and  2b  when  taken  together  dem¬ 
onstrate  that  for  properly  chosen  inert  gas-vapor  working 
fluids  the  idealized  efficiency  and  normalized  cooling  power 
will  be  increased  compared  to  a  dty  working  fluid  as  long  as 
the  normalized  temperature  gradient  lies  between  zero  and 

Having  thus  described  a  preferred  embodiment  of  the 
invention  in  sufficient  detail  to  enable  those  skilled  in  the  art 
to  make  dnd  use  the  invention,  it  will  nevertheless  be 
appreciated  that  numerous  variations  and  modifications  of 
the  illustrated  embodiment  may  be  made  without  departing 
from  the  spirit  of  the  invention,  and  it  is  intended  that  the 
invention  not  be  limited  by  the  above  description  or  accom¬ 
panying  drawings,  but  that  it  be  defined  solely  in  accordance 
with  the  appended  claims. 

We  claim: 

1.  A  thermoacoustic  refrigeration  device,  comprising: 

a  housing  arranged  to  contain  a  working  fluid; 

an  acoustic  driver  arranged  to  induce  acoustic  waves  in 
aid  working  fluid;  and 

a  thermal  stack  situated  within  the  housing  and  arranged 
to  transfer  heat  energy  from  the  working  fluid  to  a  hot 
heat  exchanger,  and  to  supply  heat  energy  from  a  cold 
heat  exchanger  to  the  working  fluid, 

wherein  the  working  fluid  is  a  gas-vapor  mixture,  and 

wherein  the  stack  is  wettable  by  condensed  vapor  such 
that  vapor  condenses  on  the  stack  during  a  refrigeration 
cycle  in  order  to  expedite  said  transfer  of  each  energy 
from  the  working  fluid  to  the  hot  heat  exchanpr. 

2.  A  thermoacoustic  refrigeration  device  as  claimed  in 
claim  1,  wherein  said  stack  is  a  finely  divided  structure 
composed  of  a  solid  material  and  having  a  length  that  is  less 
than  a  wavelength  of  said  waves. 

3.  A  thermoacoustic  refrigeration  device  as  claimed  in 
claim  1,  wherein  said  housing  includes  a  hard  termination  at 
an  end  opposite  an  end  where  the  acoustic  driver  is  situated, 
and  wherein  said  acoustic  waves  are  standing  waves. 

4.  A  thermoacoustic  refrigeration  device  as  claimed  in 
claim  1,  wherein  said  working  fluid  is  selected  to: 

a.  maximize  the  amount  of  heat  carried  by  mass  relative 
to  the  amount  of  heat  carried  by  thermal  effects,  rep¬ 
resented  by  Cp  defined  by  the  relationship: 


}  y 

CpPo  y-\ 


(w 


““  Sliquid )i 


where  p,"’  is  the  mass  density  of  the  vapor  in  the  mixture, 
is  the  mass  density  of  the  mixture,  is  the  number  density 
of  the  mixture,  n/  is  the  number  density  of  the  inert  gas,  s^£, 
is  the  entropy  per  unit  mass  of  the  mixture  ^uquid  entropy 
per  unit  mass  of  the  condensed  liquid,  and  c^  is  the  heat 
capacity  at  constant  pressure  per  unit  mass  of  the  mixture 
and, 

b.  minimize  the  heat  transfer  coefficient  9  defined  by  the 
relationship: 


where  y  is  the  ratio  of  the  specific  heats  of  the  gas  and  of  the 
mixture;  1  is  the  latent  heat  of  vaporization  of  the  mixture  per 
mole,  is  the  ambient  temperature,  and  is  the  universal 
gas  constant.  This  parameter  controls  the  critical  gradient 
10  2nd  therefore  the  stack  length  necessary  to  produce  a  given 
temperature  difference. 

5.  A  thermoacoustic  refrigeration  device  as  claimed  in 
claim  1,  further  comprising  return  means  for  transporting 
condensed  vapor  from  one  end  of  said  stack  to  the  other. 
15  6.  A  thermoacoustic  refrigeration  device  as  claimed  in 

claim  5,  wherein  said  remm  means  is  an  external  wick. 

7.  A  thermoacoustic  refrigeration  device  as  claimed  in 
claim  5,  wherein  said  return  means  is  an  external  pump. 

8.  A  thermoacoustic  refrigeration  device  as  claimed  in 
claim  1,  wherein  said  condensed  vapor  is  arranged  to  return 
to  an  opposite  side  of  the  stack  by  force  of  gravity. 

9.  A  thermoacoustic  refrigeration  method,  comprising  the 
steps  of: 

25  acoustically  driving  a  working  fluid  to  cause  waves  in  the 
working  fluid  to  transport  heat  energy  from  one  side  of 
a  thermal  stack  to  the  other  side; 

evaporating  a  vapor  in  the  working  fluid  to  increase  heat 
transfer  to  the  working  fluid  at  one  end  of  the  stack;  and 
30  condensing  the  vapor  to  increase  heat  transfer  from  the 
working  fluid  at  a  second  end  of  the  stack. 

10.  A  thermoacoustic  refrigeration  method,  comprising 
the  following  refrigeration  cycle: 

a.  an  acoustic  wave  applied  by  an  acoustic  driver  causes 
35  a  parcel  of  gas  in  a  working  fluid  to  undergo  translation 

along  a  thermal  stack  and  consequent  acoustic 
compression,  thereby  decreasing  the  parcel* s  volume 
and  increasing  its  temperature; 

b.  the  decreased  volume  and  increased  temperature 
40  increases  the  partial  pressure  of  the  vapor  within  the 

parcel; 

c.  the  parcel  then  slows,  stops,  and  reverses  its  transla¬ 
tional  motion,  while  at  the  same  time  exchanging  heat 
and  vapor  with  the  stack  as  a  result  of  the  parceFs 

45  increased  temperature  relative  to  the  stack; 

d.  at  the  time  of  reversal,  the  increased  partial  pressure 
relative  to  the  vapor  pressure  at  the  stack  wall  causes 
vapor  to  condense  from  the  parcel  to  the  adjacent  stack 
plate; 

50  e.  the  gas  parcel  then  undergoes  acoustic  rarefaction  and 
is  translated  back  past  the  ambient  position,  increasing 
its  volume  and  decreasing  its  temperature; 

f  the  acoustic  rarefaction  in  turn  causes  a  decrease  in 
partial  pressure  of  vapor  within  the  parcel; 

55  g.  the  parcel  again  slows,  stops  and  reverses  its  transla¬ 
tion^  motion  while  exchanging  heal  and  vapor  with  the 
stack,  this  time  absorbing  heal  from  the  stack;  and 

h.  since  the  partial  pressure  of  the  vapor  in  the  parcel  is 
lower  than  the  vapor  pressure  at  the  slack  wall,  the 
60  vapor  wDl  evaporate  to  the  parcel  from  the  liquid  layer 
coating  the  stack. 

♦  ♦  4:  *  + 


An  Impedance  Tube  Investigation  of  Porous  Materials 
in  the  Presence  of  a  Multi-phase  Working  Fluid 


A  Thesis 
Presented  for  the 
Master  of  Science 
Degree 

The  University  of  Mississippi 

Dan  Brown 
August  2003 


Abstract 


An  impedance  tube  is  used  to  measure  the  acoustical  impedance  of  a  porous  sample 
(thermoacoustic  stack)  with  a  working  fluid  in  multiple  phases.  The  acoustic  medium  is 
a  mixture  of  an  inert  gas,  helium,  and  a  condensable  fluid,  either  n-butane  or  acetone. 

The  results  are  compared  to  predictions  of  first  order  acoustical  theory.  Two-microphone 
impedance  measurements,  impedance  translation,  and  problems  of  fluid  handling  are 
discussed. 
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I.  Introduction 

Thermoacoustic  refrigerators  will  never  achieve  widespread  use  unless  their 
efficiency  can  be  increased.  The  inefficiency  of  these  machines  is  not  limited  by 
current  engineering,  but  is  inherent  to  this  method  of  heat  pumping.  Recent  theoretical 
developments  [Raspet  (2002),  Slaton  (2002)]  look  to  increase  both  the  efficiency  and 
heat  pumped  in  the  thermoacoustic  cycle  by  adding  a  condensable  vapor  to  the  inert 
working  gas.  This  additional  vapor  undergoes  phase  changes  between  liquid  and  gas 
during  an  acoustic  cycle.  Exploitation  of  the  latent  heat  accompanying  the  phase 
changes  may  allow  increases  in  both  efficiency  and  power  density  for  temperature 
gradients  that  provide  significant  heat  transfer  conditions. 

This  theory  requires  the  surface  of  the  stack  be  coated  with  a  thin  fluid  layer. 

We  propose  a  method  for  wetting  the  stack.  Will  the  presence  of  a  saturated 
condensable  vapor  form  a  thin  layer  over  all  the  surfaces  of  a  stack?  If  so,  does  the 
proposed  theory  correctly  describe  acoustic  propagation  in  an  inert  gas/vapor  mixture? 
How  does  the  presence  of  vapor  such  as  acetone,  isopropanol,  or  n-butane  affect  the 
elements  of  an  acoustic  system? 

To  answer  these  questions,  we  built  a  two-microphone  impedance  tube  capable 
of  measuring  stacks  in  the  presence  of  condensable  vapors.  We  directly  measured  the 
impedance  at  the  face  of  the  stack  and  compared  this  to  the  predicted  impedance  at  the 
same  location.  A  model  using  the  theory  proposed  by  Slaton,  et.  al.,  generated  the 
predicted  impedance. 
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We  begin  by  outlining  the  basic  theoiy  of  thermoacoustics  with  a  vapor  present 
and  the  equations  relevant  to  our  impedance  tube  measurements.  An  explanation  of  the 
two-microphone  method  of  measuring  acoustic  impedance  follows,  as  well  as  the 
description  of  the  system  we  have  built.  We  conclude  with  the  data  collected  and 
discuss  what  we  learned  and  ideas  for  future  work. 

II.  Theory  -  Wet  Thermoacoustics 

In  this  section,  we  establish  a  qualitative  xmderstanding  of  the  importance  and 
mechanism  of  thermoacoustics  in  the  presence  of  a  condensable  working  fluid.  After 
this  is  covered,  we  want  to  study  the  equations  governing  acoustic  propagation  in  this 
type  of  working  fluid.  These  equations  illustrate  that  the  effects  of  a  condensable 
working  fluid  will  create  a  detectable  change  in  the  measured  impedance  of  a  stack. 
Finally,  we  discuss  the  code  used  to  generate  the  theoretical  data  to  which  the  measured 
data  is  compared. 

The  details  of  wet  thermoacoustics  can  be  daunting;  a  simplified  explanation 
will  illuminate  the  important  concepts.  In  the  Lagrangian  picture  of  thermoacoustic 
refrigerators  [Swift  (1988)],  Fig.  (1),  we  have  a  parcel  of  gas  near  some  surface  whose 
position  and  pressure  oscillate.  At  each  extreme  of  its  displacement  the  parcel 
-experiences  either  a  maximum  or  minimum  acoustic  pressure.  The  compression  and 
rarefaction  cause  a  temperature  difference  between  the  gas  parcel  and  surface.  Heat 
will  be  absorbed  from  the  surface  by  the  rarefied  parcel  and  deposited  to  the  surface  by 
the  compressed  parcel.  This  develops  a  temperature  difference  along  the  surface  near 
the  oscillating  parcel.  This  is  shown  in  Figure  (1)  as  heat  is  pumped  from  left  to  right. 
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stack  Wall 

Figure  (1):  A  Lagrangian  model  for  a  dry  thermoacoustic  refrigerator.  Compression 
and  rarefaction  of  a  gas  parcel  in  the  acoustic  cycle  allow  for  heat  pumping  due  to 
changes  in  gas  temperature. 
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Although  the  practical  design  of  a  thermoacoustic  refrigerator  is  simple, 
applications  are  hampered  by  their  inherent  inefficiencies.  The  addition  of  a 
condensable  vapor  to  these  machines  may  increase  their  efficiency  as  well  as  power 
density.  For  temperature  gradients  near  the  critical  temperature  gradient, 
thermoacoustic  refrigerators  are  extremely  efficient  (Figure  2),  but  these  operating 
conditions  are  not  useful  because  of  the  small  power  density  associated  with  efficient 
heat  pumping.  For  significant  heat  pumping,  as  required  by  realistic  machines,  the 
efficiency  is  lower,  and  the  temperature  gradient  may  be  less  than  50%  of  the  critical 
temperature  gradient.  Looking  at  the  efficiency  and  power  density  of  a  helium 
refrigerator  at  F  =  0.5 ,  we  see  an  equal  or  better  efficiency  can  be  attained  using  a 
helium/n-butane  mixture  over  0.25  <  F  <  0.40 ,  and  the  power  density  is  superior  to  that 
of  helium  over  the  entire  range. 

A  compromise  between  efficiency  and  cooling  power  must  be  reached  in  the 
design  of  thermoacoustic  refrigerators.  The  addition  of  condensable  working  fluids 
creates  regions  over  certain  temperature  gradients  where  increases  in  both  efficiency 
and  cooling  power  are  possible. 
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Efficiency 


r  =  (GRAD  T)  /  (GRAD  T)cRn 


Figure  (2):  Comparison  of  normalized  coefficient  of performance  and  normalized 
cooling  power  for  wet  and  dry  thermoacoustic  refrigerators.  For  a  given  temperature 
gradient,  it  is  possible  to  increase  both  cooling  power  and  efficiency  over  that  of  a  dry 
system  [Slaton  (2002)]. 
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Figure  (3):  A  Lagrangiart  model  for  a  wet  thermoacousHc  refrigerator.  Mass  diffusion 
is  present  along  with  sensible  heat  pumping.  Evaporation  and  condensation  at  the 
liquid  layer  increases  the  heat  pumped  with  each  acoustic  cycle. 

In  the  previous  example  of  dry  thermoacoustics,  the  only  working  fluid  present 
in  the  system  is  an  inert  gas  or  a  mixture  of  inert  gasses.  In  a  wet  thermoacoustic 
refrigerator,  the  working  fluid  consists  of  a  mixture  of  inert  gas  and  condensable  vapor. 
The  vapor  is  added  to  the  inert  working  gas  until  the  system  is  saturated.  This  occurs 
when  the  rate  of  liquid  evaporation  equals  the  rate  of  vapor  condensation,  or  when  the 
vapor’s  partial  pressure  equals  the  liquid’s  vapor  pressure.  It  is  assumed  that  once 
saturation  is  reached,  all  surfaces  have  a  thin  layer  of  liquid  present  on  them.  The 
presence  of  both  the  liquid  and  gas  phase  in  equilibrium  is  the  key  to  wet 
thermoacoustics. 

During  the  acoustic  cycle,  the  equilibrium  is  disturbed.  The  vapor’s  partial 
pressure  will  depart  from  the  liquid’s  vapor  pressure,  resulting  in  evaporation  or 
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condensation.  In  a  Lagrangian  picture  similar  to  that  for  the  dry  case  (Figure  3),  a 
concentration  gradient  develops  between  the  rarefied  and  compressed  parcels  and  the 
surrounding  vapor.  The  result  is  mass  diffusion.  As  the  gas  parcel  is  rarefied,  the 
partial  pressure  of  the  vapor  in  that  parcel  is  reduced  below  the  partial  pressure  of  the 
vapor  above  the  adjoining  liquid  layer.  A  concentration  gradient  develops  and  vapor 
will  pass  into  the  parcel.  The  vapor  that  passes  into  the  parcel  must  be  replaced  in  order 
to  maintain  equilibrium  with  the  liquid  layer.  As  a  result,  some  of  the  liquid  evaporates. 
The  surface  now  loses  the  latent  heat  of  vaporization  in  addition  to  the  sensible  heat  lost 
due  to  adiabatic  expansion  of  the  gas  parcel.  Likewise,  a  compressed  parcel  develops  a 
concentration  gradient  because  the  partial  pressure  of  the  parcel’s  vapor  is  greater  than 
that  of  the  surrounding  vapor.  Vapor  diffuses  out  of  the  parcel.  Since  the  vapor  at  the 
surface  is  in  equilibrium  with  the  liquid,  the  diffused  vapor  must  condense. 
Condensation  causes  the  surface  to  absorb  the  latent  heat.  Hence,  we  get  additional  heat 
pumping  for  an  acoustic  cycle.  Firstly,  we  have  the  adiabatic  heat  pumping  due  to  the 
standard  model  of  thermoacoustics,  but  we  also  add  the  heat  pumping  caused  by  mass 
diffusion. 

Now  that  we  have  a  conceptual  outline  of  wet  thermoacoustics,  we  will  look  at 
the  mathematics  of  acoustic  propagation  in  saturated  vapors.  We  are  looking  to  see 
how  the  presence  of  vapor  and  accompanying  mass  diffusion  will  affect  the  impedance 
measured  in  our  system.  Looking  at  the  equations  that  govern  the  acoustics  of  dry  and 
wet  gases,  we  see  many  similarities.  The  form  of  the  wave  equation  is  the  same  for 
both  wet  and  dry  pores  [Raspet  (2002)]. 
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(1) 


d^P  ^  dP  ,2_  _ 

dz^  dz 

Where  for  dry  pores  the  attenuation  factor,  not  to  be  confused  with  the  imaginary  part  of 
k  which  is  the  attenuation,  and  wave  number  are  given  by 
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And  for  wet  pores. 
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where, 

Y-\  I 


<^  =  - 


(6) 


r 

is  the  latent  heat  parameter.  For  the  square  pores  being  investigated  we  have  [Amott 
(1991)] 


FiX)  =  ^  ' 


(7) 


where. 


4.W  =  l+ 


(8) 
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Reducing  equation  (1)  to  a  pair  of  first  order  differential  equations  illuminates  the  way 
in  which  the  pressure  and  volume  velocity  evolve  in  the  system  [Raspet  (2002)]. 
dPAc  _  i(op  .Q. 

€I^^^F(X)k'‘P^-2aU^  (10) 

fife  lODp 

The  imaginary  part  of  k  is  the  loss  term  due  to  both  heat  and  mass  transport, 
and  a  represents  a  thermoacoustic  gain  or  attenuation  factor.  This  set  of  equations  is 
sufficient  to  model  the  acoustic  behavior  for  a  stack  if  the  geometry  of  the  system  is 
known. 

For  this  study,  the  resonator  and  stack  are  at  constant  temperature  and  the 
pressure  amplitudes  used  are  not  sufficient  to  create  a  temperature  gradient; 
consequently,  or  =  0 .  Therefore,  the  effect  of  the  vapor  will  only  be  detectable  in 
changes  to  the  complex  wave  number. 
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Figure  (4):  The  influence  of  mass  diffusion  on  the  wavenumber  in  the  stack  and  the 
impedance  in  a  model  of  our  system.  The  impedance  shown  is  a  theoretical  prediction 
of  the  impedance  we  expect  to  model  The  lines  denoted  “with  mass  diffusion  ”  use 
equation  (5)  to  calculate  the  wavenumber  inside  the  stack.  The  lines  denoted  “without 
mass  diffusion  ”  set  mass  diffusion  coefficient,  to  zero  in  the  calculation  of  the  gas 
properties.  The  effect  of  setting  Z),2  =  0  iy  to  calculate  the  wavenumber  based  on  the 
gas  properties  of  the  helium  +  n-butane  mixture  without  including  the  effects  of  mass 
transport.  That  is,  we  model  a  theoretical  gas  that  has  the  same  thermophysical 
properties  as  helium  +  n-butane  but  does  not  change  phase  during  the  acoustic  cycle  to 
obtain  a  “without  mass  diffusion  ”  curve. 

In  Figure  (4),  we  see  a  considerable  change  in  the  imaginary  part  of  the 
wavenumber  corresponding  to  an  increase  loss  during  the  acoustic  cycle.  The  real  part 
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of  the  wavenumber  does  not  change  much,  so  we  expect  the  resonance  frequency  to 
remain  almost  constant.  The  impedance  graphs  in  Figure  (4)  indicate  that  mass 
diffusion  significantly  changes  the  amplitude  of  the  measured  impedance  at  the  face  of 
the  stack.  As  suspected,  the  resonance  frequency  will  not  change  dramatically.  In  order 
to  identify  this  shift,  our  system  must  be  able  to  resolve  changes  of  1  Hz. 

The  theory  to  which  we  compared  our  dry  measurement  was  coded  using  work 
published  by  Amott  [Amott  (1991)].  We  used  the  known  value  of  the  impedance  of  a 
hard  end  with  thermal  dissipation  [Pierce  (1991)]  and  the  intrinsic  impedances  of  the 
tube  and  stack  to  predict  the  impedance  at  the  face  of  the  stack  using  the  impedance 
translation  theorem  [Pierce  (1991)].  See  Appendix  B  for  these  equations.  This  was 
compared  to  the  impedance  measured  at  flie  face  of  the  stack  (Figures  5  &  6). 


Excite  Acoustic  Amplify  Signals  Obtain  Transfer  Calculate  Impedance 

Waves  In  - ►  From  Kezoresistive  - ►  Function  Brom  - ►  At  Stack  face  From 

Impedance  Tube  TransduCeis  Signal  Analyzer  Transfer  Rjiiction 


Figure  (5):  Flow  chart  for  acquisition  of  data  from  the  experimental  system. 
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Calculate  Impedance  Translate  Hard  End  Translate  Impedance  Compare  To 

Of  Hard  End  With  - ►  Impedance  Through - Throu^  Stack  - •-  Measured  Impedance 

Thermal  Dissipation  Pipe  Using  ITT  Using  ITT  At  Stack  Face 


Figure  (6):  Flow  chart  for  theoretical  prediction  of  impedance  at  stack face  using  the 
impedance  translation  theorem  (ITT). 

The  wet  theory  code  was  produced  in  the  same  manner.  Of  primary  importance 
in  the  calculation  of  the  wet  impedance  are  the  gas  properties,  which  were  obtained 
from  gas  dynamic  theory  [Hirschfelder  (1964)].  Dr.  William  Slaton  provided  all 
theoretical  calculations  for  the  cases  of  wet  stacks. 
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Acetone  2-Propanol  n-Butane 


Temperature  (K) 

3.000E+02 

3.000E+02 

3.000E+02 

li  i  ll  Mill  tixi  :4!/ii 

2.321E+00 

1.764E+00 

6.684E+00 

cp  mixture  (J/kg  K) 

3.937E+03 

4.877E+03 

2.281  E+03 

Thermal  Conductivity  mixture  (W/m  K) 

1.379E-01 

1.539E-01 

5.958E-02 

sound  speed  mixture  (m/s) 

8.222E+02 

9.596E+02 

4.433E+02 

Viscosity  mixture  (pPa  s) 

1.975E+01 

2.009E+01 

1.364E+01 

mole  fraction  gas 

9.936E-01 

7.689E-01 

mole  fraction  vapor 

6.352E-03 

2.31  IE-01 

mole  density  gas 

3.920E+02 

4.024E+02 

3.114E+02 

mole  density  vapor 

1.297E+01 

2.572E+00 

9.357E+01 

Ratio  of  specific  heats  for  mixture 

1.554E+00 

1.608E+00 

1.301  E+00 

D12  (m^/s) 

4.538E-06 

4.447E-06 

3.515E-06 

Prandtl  Number 

6.366E-01 

5.221  E-01 

Schmidt  Number 

1.875E+00 

2.561  E+00 

5.804E-01 

Phi  -  Latent  heat  parameter 

MiSiHfB! 

6.687E+00 

1.930E+00 

Thermal  Diffusion  Ratio 

-6.485E-03 

-1.273E-01 

Table  (1 ):  Properties  of  Helium  +  vapor  at  10  bar. 

ni.  Theory  -  Measurement 

In  1980,  Chung  and  Blaser  [Chung  &  Blaser  (1980)]  described  a  technique  of 
measuring  impedance  by  using  broadband  noise  and  two  transducers  mounted  in  an 
impedance  tube.  Through  a  comparison  of  the  two  measured  pressures,  the  acoustic 
impedance  is  determined  at  any  other  position  in  the  system. 

Comparing  the  two  signals,  we  take  the  ratio  of  the  cross-spectrum  5,2  (/)  to 
the  auto-spectrum  5„  (/) .  This  ratio  is  also  known  as  the  transfer  function  (/)  •  It 
is  from  the  transfer  function  the  impedance  is  calculated. 


HM)- 


Snif) 

Sn(f) 


(11) 
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(12) 


Z  _  .  sin[A^(/-5)]-//i2  sin(A:/) 
pc  cos[k(l  -  ^)]  -  i/,2  cos{kl) 

Where  s  is  the  separation  of  the  transducers,  k  is  the  wavenumber,  and  /  is  the  distance 
from  the  first  transducer  to  the  location  where  the  impedance  is  measured.  In  this  work 
this  is  taken  at  the  face  of  the  stack  towards  the  driver.  We  use  the  real  wavenumber 

k-  —  neglecting  the  small  attenuation,  a  =  0.0021  Np/m,  between  the  microphone 
c 

and  the  sample. 

The  two-transducer  technique  has  several  advantages  over  other  methods  of 
measuring  acoustic  impedance.  Modem  digital  signal  analyzers  are  capable  of  directly 
measuring  the  transfer  function,  thereby  greatly  reducing  the  need  for  complex  signal 
processing  routines  to  be  coded.  Additionally,  the  broadband  nature  of  this 
measurement  allows  for  rapid  data  acquisition  over  the  frequency  range  of  interest. 
However,  measurement  of  the  transfer  function  can  be  corrupted  by  error  and 
nonlinearity  in  the  pressure  measurement.  Detail  of  the  measurement  system  and 
calibration  procedure  can  be  found  in  Appendix  A. 
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IV 


Experimental  Setup 


Figure  (7):  Schematic  of  System 


Figure  (8):  Photo  of  system 


20 


Figure  (9):  Sample  holder  and  transducers 

a.  General  Description 

Figure  (7)  is  a  schematic  of  the  impedance  tube  used  in  this  work.  The  inner 
diameter  of  the  tube  is  0.1016  m.  The  total  length  is  about  2.5  m.  In  figure  (8),  the 
mounting  r^od  for  the  speaker  can  be  seen  protruding  from  the  leftmost  flange.  The 
speaker  used  is  a  0.1016  m  driver  clearanced  to  fit  in  the  tube  and  powered  by  an 
external  amplifier.  Mounted  on  the  right  end  of  the  tube  is  a  blind  flange  used  as  the 
rigid  termination.  The  tube  and  all  flanges  are  made  of  stainless  steel.  The  sample  is 
housed  in  the  short  section  shown  in  Figure  (9).  Just  to  the  left  of  this  section,  the  two 
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transducers  are  visible.  Fine  gauge  stainless  tubing  is  connected  to  the  vent  port  of  the 
transducers  in  order  to  maintain  the  reference  side  of  the  transducers  at  the  static 
pressure  of  the  tube.  Once  the  tube  is  at  operating  pressure,  a  valve  isolates  the 
reference  side  of  the  transducer  from  the  sound  field.  This  allows  the  transducers  to  see 
a  varying  pressure  only  on  the  face  exposed  to  the  resonator.  To  the  right  of  the  sample 
holder,  we  have  mounted  a  high  precision  thermometer  as  well  as  the  fittings  necessary 
to  fill  the  tube  with  gas  and  vapor, 
b.  Corrosive  Vapors 

As  previously  stated,  wet  thermoacoustics  uses  a  mixture  of  inert  gas  and  vapor 
as  the  working  fluid.  This  vapor  must  be  present  in  quantities  sufficient  to  coat  all 
surfaces  within  the  resonator  with  a  thin  layer  of  liquid.  Typical  candidates  for  the 
condensable  component  are  hydrocarbons,  alcohols,  freon,  and  water. 

One  of  the  concerns  when  starting  this  project  was  the  durability  of  the  materials 
to  various  liquids  and  vapors.  Acetone  and  n-butane  are  effective  at  dissolving  bonding 
agents  as  well  as  some  seals.  This  presents  a  problem  with  PTFE  seals,  silicone 
lubricants,  and  bonding  agents.  The  glues  used  in  the  assembly  of  our  drivers  were 
particularly  susceptible  to  exposure  to  acetone  and  butane.  Fig.  (10).  Working  fairly 
.  quickly,  in  less  than  a  few  hours,  allowed  collection  of  useful  data  before  the 
disintegration  of  the  driver. 
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Figure  (10):  Damaged  driver  (bonding  agents  dissolved)  that  was  exposed  to  n-butane. 


This  driver  was  usable  for  approximately  150  minutes  in  a  helium/butane  environment, 
an  obviously  unacceptable  result  for  any  system  designed  for  continuous  use. 
c.  Transducer  separation 

We  found  that  multiple  transducer  separations  are  necessary  for  the  range  of 
wavelengths  we  would  encounter.  The  maximum  separation  s  is  given  by 


(12) 
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where  is  the  minimum  frequency  to  be  investigated.  This  simply  limits  the 

separation  to  ensure  that  you  have  at  most  a  half  wavelength  between  microphones 
since  the  reduction  algorithm  assumes  zero  or  one  pressure  node  between  transducers. 
A  minimum  separation  is  not  as  clearly  defined.  Narrowly  separated  transducers  do  not 
resolve  differences  in  phase  and  amplitude  well,  especially  when  dealing  with  long 
wavelengths.  We  must  strike  a  balance  for  the  separation  to  be  used.  For  any  given 
mixture,  we  want  maximum  separation  for  amplitude  and  phase  resolution,  but  we 
cannot  exceed  the  limits  placed  by  Eq.  12. 

The  mixtures  we  investigated  have  speeds  of  sound  that  vary  by  a  factor  of 
three,  and  we  investigate  them  over  a  fi-equency  range  that  spans  an  order  of  magnitude. 
With  wavelengths  ranging  from  0.3  m  to  10.0  m,  multiple  transducer  separations  were 
necessary.  We  settled  on  three  separations:  0.0381m,  0.1143  m,  and  0.1524  m.  This 
range  permits  flexibility  for  gas/vapor  selection, 
d.  Transducer  Selection 

Two  factors  were  considered  in  the  transducer  selection:  sensitivity  and  size. 
Initial  measurements  were  made  that  showed  pressure  amplitudes  on  the  order  of  1 .0 
Pa.  With  acoustic  signals  5  orders  of  magnitude  smaller  than  anticipated  maximum 
static  pressures,  differential  transducers  are  necessary.  Both  sides  of  a  differential 
transducer’s  membrane  are  exposed  to  the  static  pressure,  allowing  for  a  very  high 
sensitivity  since  the  membrane  only  reacts  to  acoustic  pressures. 

The  size  of  the  transducer  determines  how  it  will  be  mounted  in  the  tube  wall. 

In  this  system,  the  tube  wall  is  0.25”  thick,  and  the  pipe  has  an  outer  diameter  of  4.5”. 
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The  curvature  and  thickness  of  the  pipe  dictated  that  a  narrow  transducer  be  used  if 
extraordinary  means  of  mounting  are  to  be  avoided.  An  Endevco  model  8510-Bl 
satisfies  the  above  parameters.  It  has  a  differential  pressure  range  of  1  psi  and  a  full- 
scale  output  of  ~200  mV  with  10  V  excitation  potential.  In  addition  to  satisfying  the 
above  characteristics,  they  are  relatively  inexpensive,  so  that  damage  due  to  over 
pressurization  is  not  financially  prohibitive. 


Figure  (11):  Endevco  Model  8510B-1  piezoresistive pressure  transducer 


25 


1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - ] - 1 - 1 - 1 - 1 - ] - 1 - 1 - 1 - 1 - 1 - 1 - 1 

-1.0  -0.5  0.0  0.5  1.0 

Pressure  (PSIG) 

Figure  (12):  Examination  of  the  Endevco  8510B-1  response  to  DC  pressure 
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e.  Data  Sets 

Data  extracted  from  the  spectrum  analyzer  consists  of  the  FFT  of  each 
transducer’s  signal,  the  real  and  imaginary  parts  of  the  transfer  function  between 
transducers,  and  the  coherence  of  the  signal.  The  coherence  is  helpful  in  determining  if 
the  signal-to-noise  ratio  is  adequate.  Figure  (1 1)  shows  coherence  indicative  of  a  poor 
signal-to-noise  ratio. 


Figure  (13):  Coherence  when  proper  signal-to-noise  ratio  is  not  obtained.  This 
particular  plot  is  of  helium  at  0.10  bar  of pressure  at  room  temperature. 
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Pressure  Amplitude  (Volts) 


The  FFT  of  the  two  channels  for  this  case  is  shown  in  Figure  (14). 


Figure  (14):  FFT  of  data  set  with  insufficient  signal-to-noise  ratio. 


28 


Coherence 


Data  shown  in  Figures  (15),  (16)  demonstrates  sufficient  signal-to-noise  ratios; 


consequently,  the  coherence  is  close  to  1. 


Figure  (15):  Example  of  proper  coherence.  Data  set  is  helium  at  10. 13  bar  with  a  600 
PPf  stack  at -295  K. 
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Presssure  Amplitude  (Volts) 


Figure  (16):  FFT  of  data  set  with  sufficient  signal-to-noise  ratio. 
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V.  Data  &  Analysis 


Fig.  (1 7):  Celcor  stack.  This  is  the  same  material  as  the  stack  we  investigated.  This 
stack  has  600 pores  per  square  inch  and  is  made  of  ceramic. 


31 


We  measured  the  specific  acoustic  impedance  of  Celcor  stacks  in  both  inert  gas 
and  inert  gas/yapor  mixtures  -  see  Figure  (17).  The  stack  investigated  has  600  pores  per 
square  inch,  and  a  porosity  of  78%.  All  data  sets  were  taken  at  temperatures  in  the 
range  of 294  K  -  297  K.  For  the  data  using  only  inert  gas,  we  investigated  helium  and 
argon  at  1  bar  and  10  bars.  All  of  the  following  plots  are  of  the  specific  acoustic 
impedance  at  the  face  of  the  stack. 
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Imaginary  Imedance  (Z/pc)  Real  Imedance  (Z/pc) 
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Imaginary  Imedance  (Z/pc)  Real  Imedance  {ZIpc) 
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Figure  (18):  The  above  graphs  compare  the  measured  and  calculated  impedance  at  the 
face  of  the  stack  for  helium  and  argon  at  1.01  bar  and  10.13  bar. 

In  all  of  these  sets,  we  see  that  the  data  and  theory  have  the  same  form;  however, 
there  is  a  discrepancy  in  both  jfrequency  and  amplitude.  Remembering  the  condensable 
vapor  causes  small  shifts  in  the  resonance  frequency,  this  discrepancy  is  unacceptable. 
We  must  either  eliminate  the  cause  of  this  dissimilarity  or  characterize  it.  Examination 
of  the  frequency  difference  in  the  argon  @10  bar  graph  shows  characterization  is 
difficult.  The  shift  is  not  uniform  or  linear.  At  the  first  resonance  peak,  a  difference  of 
approximately  5  Hz  exists  between  data  and  theory,  while  subsequent  peaks  show  shifts 
of  10  Hz.  Similar  problems  present  themselves  in  attempts  to  characterize  the 
amplitude  discrepancies. 

The  source  of  the  differences  went  undiscovered.  Three  possibilities  exist:  the 
data  taken  does  not  represent  the  impedance  at  the  face  of  the  stack;  the  theory  used 
does  not  predict  the  intrinsic  impedance  of  the  stack  correctly;  or  our  assumptions 
regarding  the  resonator  are  not  correct. 

First,  we  examined  the  experimental  measurement  of  the  impedance  at  the  face 
of  the  stack.  All  the  appropriate  dimensions  were  remeasured,  and  the 
amplifier/transducer  system  was  recalibrated.  The  recalibration  did  not  produce 
significant  changes  in  the  intrinsic  transfer  function,  nor  did  the  remeasurement  reveal 
any  errors. 

As  an  alternative  to  using  the  intrinsic  transfer  function  calibration  (Appendix 
A),  it  is  possible  to  employ  sensor  switching  to  eliminate  gain  and  phase  errors 
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introduced  by  the  amplifier/transducer  system.  We  measured  the  impedance  using  this 
technique  and  compared  it  to  data  taken  using  the  calibration  technique.  The  agreements 
between  these  two  data  sets  led  us  to  believe  that  the  experimental  measurements 
accurately  represent  the  impedance  at  the  face  of  the  stack. 

Next,  we  addressed  the  theory  and  its  application.  The  theory  published  by 
Amott  was  coded  in  C-H-.  The  output  of  this  code  was  compared  to  an  identical  model 
prepared  in  deltaE™.  Good  agreement  existed  between  these  two  simulations  as  shown 
in  Fig.  (19). 
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Figure  (19):  Comparing  our  code  to  deltaE  to  confirm  it  correctly  predicts  the 
impedance  for  the  programmed  geometry. 


Now  we  turned  to  our  assumptions  in  modeling  sections  of  the  impedance  tube 
other  than  the  stack.  The  model  coded  is  that  of  a  hard  end  followed  by  a  duct,  then  a 
stack.  We  integrate  the  pressure  and  volume  velocity  equations  from  the  termination 
through  to  the  face  of  the  stack  and  this  gives  a  measure  of  the  impedance.  We  do  take 
into  consideration  the  thermal  dissipation  at  the  hard  end  and  the  attenuation  at  the  tube 
wall.  We  do  not  take  into  account  any  compliance  that  may  come  about  due  to  the  gas 
inlet  tubes,  the  transducer  reference  pressure  plumbing,  or  the  gaskets  between  the 
flanges.  These  are  all  things  that  could  be  the  cause  of  die  problems  we  are  facing. 
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While  it  may  be  possible  to  model  these  additional  compliances,  determination  of  their 
true  influence  may  be  difficult  to  ascertain. 

Another  source  of  error  may  be  in  the  modeling  of  the  stack  material  itself.  The 
walls  of  the  ceramic  stack  we  are  using  are  porous  and  this  may  affect  the  complex 
wavenumber  within  the  stack.  Amott  found  wall  porosities  as  high  as  49%  are  needed 
to  model  Celcor  properly  [Amott  (1991)].  We  conducted  an  experiment  to  determine 
both  the  wall  porosity  and  the  porosity  due  to  the  main  pores.  We  immersed  the  stack 
in  water,  removed  it,  and  used  compressed  air  to  blow  the  water  out  of  the  main  pores. 
The  mass  of  the  stack  changed  by  86  g  once  the  walls  were  wet.  Assuming  all  86  g  are 
contained  in  the  walls,  the  walls  have  86  cc  of  pores.  We  reinserted  the  wet  stack  and 
measured  how  much  water  it  displaced.  With  the  walls  wet,  the  volume  of  displaced 
water  should  equal  the  volume  of  the  walls  if  they  were  solid.  Knowing  the  total 
volume  of  the  stack,  we  calculated  a  main  pore  porosity  of  78%  and  a  wall  porosity  of 
39%.  From  Amott,  we  have  the  dispersion  relation  for  the  wavenumber  for  porous 
walls 


porous 


=  ±k 


nonporous. 


(13) 


where, 
^«1  + 


and 


(14) 


Qr  =^^(1  +  2 


R  ^ 


(15) 
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for  the  total  porosity,  is  the  average  depth  of  the  wall  pores  and,  assuming  it  is  50 


pm,  good  agreement  exists  between  data  and  theory. 
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Figure  (20):  Examining  the  influence  of  changing  the  complex  wavenumber  in  the  stack 
to  account  for  the  porous  walls  of  a  ceramic  stack. 


We  should  mention  that  the  value  used  for  ,  one  third  of  wall  thickness,  was 
chosen  to  give  good  agreement  between  data  and  theory.  Verification  of  this  number  is 
beyond  the  scope  of  this  work. 

Though  perplexed  by  the  shifts  in  data,  we  still  investigated  the  inert  gas/vapor 
mixture.  What  follows  is  one  set  of  that  data: 
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Figure  (21):  The  above  data  compares  the  theoretical  predictions  by  Slaton  to  the 
measured  data.  At  the  time  of  writing,  the  only  simulation  available  for  comparison 
was  that  of 300  Kat  10. 13  bar.  The  temperature  difference  accounts  for  a  1%  change 
in  speed  of  sound.  The  pressure  difference  will  cause  a  change  in  the  molar 
concentration  of  the  vapor  and  hence  all  the  gas  properties.  The  effects  of  this  pressure 
difference  are  not  apparent  due  to  the  complex  nature  of  the  equations. 

The  data  shown  in  Figure  (21)  show  good  agreement  in  amplitude,  and  this  is 
promising.  However,  the  errors  in  frequency  are  quite  large.  In  addition  to  this,  the 
errors  seem  to  be  proportional  to  the  frequency,  unlike  the  shift  in  the  dry  data.  How  the 
presence  of  a  fluid  layer  would  affect  the  porous  walls  of  the  stack  is  unclear,  thus  this 
was  not  investigated. 

VI.  Discussion 

The  original  goals  of  this  project  were  to  determine  if  a  stack  exposed  to  a 
saturated  vapor  develops  a  thin  layer  of  fluid  over  all  its  surfaces,  to  confirm  the  validity 
of  the  predicted  complex  wavenumber  in  the  presence  of  a  saturated  vapor,  and  to 
characterize  the  effects  vapors  would  have  on  the  elements  of  an  acoustic  system. 

The  agreement  in  amplitude  in  our  wet  data  is  hopeful,  but  the  offsets  in 
frequency  must  be  resolved  before  we  can  certify  the  wetting  of  the  stack  and  the 
changes  to  the  wavenumber.  In  order  to  accomplish  this,  several  changes  are  ^ 

recommended  for  fiitiue  work. 

The  system  should  be  redesigned  to  eliminate  as  many  potential  sources  of 
compliance  as  possible.  A  single  tube  that  is  terminated  with  a  hard  end  would  be  best. 
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Using  a  single  tube  instead  of  multiple  sections  reduces  the  possibility  of  misalignment 
as  well  as  removes  the  gasket  materials.  The  differential  microphones  should  be 
abandoned  as  well.  This  would  eliminate  the  plumbing  necessary  to  create  the 
reference  side  pressure  they  require.  Of  course,  this  presents  the  problem  of  finding  a 
transducer  that  is  sensitive  to  pressure  fluctuations  on  the  order  of  a  1  Pa,  but  is  not 
damaged  by  pressures  as  high  as  10®  Pa.  Additionally,  the  location  of  gas  and  vapor 
filling  should  be  moved  to  the  speaker  side  of  the  transducers  so  as  not  to  influence  the 
measurement. 

Elimination  of  the  wall  porosity  of  the  stack  would  be  useful  as  well,  either  by 
applications  of  varnish  [Bernard  (1996)],  or  by  selecting  a  material  that  does  not  have 
porous  walls,  such  as  stainless  steel.  This  would  allow  verification,  through  dry 
measurements,  that  the  system  correctly  measures  predicted  impedance.  However,  the 
porous  walls  may  be  useful  as  an  indicator  of  wetness  in  future  measurements.  As  the 
walls  become  coated  with  fluid  the  pores  begin  to  fill.  This  changes  the  value  of  </«-  in 
Eq.  (15),  and  the  effects  of  the  wall  porosity  on  the  complex  wavenumber.  Eventually 
the  pore  will  be  filled,  dw  =  0,  and  the  wall  porosity  no  longer  influences  the  measured 
impedance.  As  vapor  is  added  to  the  system,  and  dw  grows  smaller,  this  shift  in 
-  impedance  should  be  detectable. 

Handling  the  corrosive  vapors  presented  some  problems.  Acetone  and  n-Butane 
had  pronounced  effects  on  our  system,  the  destraction  of  the  speakers  being  the  most 
fi-equent  problem.  Finding  a  method  of  exciting  sound  waves  that  is  not  susceptible  to 
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the  vapors  being  investigated  will  be  important  if  multiple  data  runs  are  to  be  conducted 
with  each  filling  of  the  system. 

Future  work  in  this  area  is  promising.  With  the  proposed  changes  above,  I  am 
confident  that  definitive  answers  can  be  given  to  the  presented  questions. 
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Appendix  A 
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The  purpose  of  measuring  the  transfer  function  between  two  transducers  is  to 
gain  information  about  the  signals  present.  Unfortunately,  the  measurement  system 
(transducer  and  amplifier)  itself  affects  the  transfer  function.  In  addition  to  the  transfer 
function  due  to  the  signal,  ,  the  measmed  transfer  function  includes  the 

intrinsic  transfer  fimction  of  the  data  acquisition  system  . 


Measured  _  Signal  ^  ^  Intrinsic  j  ^ 

Fortunately,  accounting  for  the  nonlinear  properties  of  the  data  acquisition 
channels  is  straightforward.  To  do  this  both  channels  must  be  exposed  to  an  identical 
sound  field.  This  reduces  equation  (A  1)  to 


^Measured  ^  J  jj Intrinsic 


(A2) 


To  achieve  an  identical  sound  field,  the  transducers  were  mounted  in  radial 
symmetry  into  a  flange  (See  Picture).  This  flange  is  mounted  onto  a  pipe  and  exposed 
to  broadband  noise.  Assuming  we  have  plane  waves,  the  signal  at  each  transducer  is 
identical.  Therefore,  any  non-unity  transfer  function  will  be  due  to  the  measurement 
system  alone,  .  To  account  for  this  in  further  measurements,  is  simply 

divided  by 

It  is  worthwhile  to  note  that  i/,2  *  .  In  addition,  we  must  be  consistent  in 

determining  which  of  the  two  transducers  will  be  labeled  as  #1  or  “reference.” 

Equation  (12)  requires  the  reference  transducer  be  placed  farthest  firom  the  stack. 
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Figure  (Al):  Transducers  mounted  in  radial  symmetry  for  calibration. 


The  transfer  functions  for  our  system  follow.  The  figures  also  show  a  cubic  fit 
to  the  intrinsic  transfer  fimction,  which  was  used  to  speed  data  reduction. 
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Imaginary  Transfer  Function  (x10  )  Real  Transfer  Function 


Appendix  B 
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The  following  equations  are  used  in  conjunction  with  the  impedance  translation 
theorem  to  model  the  impedance  at  the  face  of  the  stack  with  out  the  presence  of  a 
condensable  vapor.  We  begin  with  the  impedance  of  normal  incidence  plane  wave 
reflection  at  a  solid  surface  due  to  thermal  conduction  and  viscosity  at  that  surface 
[Pierce  (1991)]. 


Z 


lit 


(Bl) 


where, 

'?,(■»)  =  (B2) 

This  is  translated  through  the  intrinsic  impedance  of  our  0.1016  m  diameter 
pipe.  This  was  calculated  using  the  boundary  layer  approximation  [Amott  (1991)]. 


QF{X)k 


where, 

F(I) 


(B3) 


(B4) 


and, 

F(A)  =  l-~(l  +  i)-  (B5) 

A> 

Now  this  impedance  can  be  translated  through  the  intrinsic  impedance  of  the 
porous  sample;  thus,  predicting  the  impedance  at  the  face  of  the  stack.  The  intrinsic 
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impedance  of  the  stack  was  calculated  using  Eqs.  B3  and  B4,  For  the  dissipation 
function  we  have 


F(;i)  = 


64^  1 


(B6) 


where, 


i;,w=i+ 


in^{m^  +n^) 

4^ 


(B7) 
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